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of  track  prediction  has  been  slow  to  oome  and  the  need  for  further 

\  . 

research  on  statistical-synoptic  approaches  exists.  This  newly 
developed  method  of  forecasting  one-  to  three-day  oyolone  notion  is  an 
improvement  over  existing  forecast  schemes. for  the  following  reasons: 

(1)  cyclones  are  stratified  based  on  their  positions  relative  to  the  500 
mb  subtropical  ridge  to  better  define  the  environmental  Influences  on 
the  cyclones;  (2)  the  72-hour  forecast  track  is  segmented  into  three 
24-hour  time  steps  to  permit  the  application  of  updated  persistence  and 
synoptic  data  reltlve  tr  each  new  cyclone  position  as  the  24-hour 
displacements  are  stepped  forward  to  the  desired  forecast  projection; 
and  (3)  the  prognostic  synoptic  data  are  Introduced  only  at  their  valid 
times  and  relative  to  the  cyclone  position  at  the  valid  time  such  that 
they  are  treated  like  analysis  data. 

A  screening  —  multiple  regression  technique. is  used  with  the 
"perfect-prog"  concept  to  Introduce  numerically  forecast  synoptic  data 
in  the  development  of  the  regression  coefficients.  Actual  forecast 


fields  are  used  to  test  the  skill  of  the  new  model  to  be  referred  to  as 


Comparison  with  operational  forecast  techniques  on  a  nearly 

homogeneous  sample  using  a  persistence-climatology  model  as  a  normalizer 

to  Indicate  relative  skill  of  the  models  as  compared  to  a  no-skill 

forecast  indioate  that  CSU84  is  six  times  as  skillful  as  the  best 

operational  model  in  the  Atlantic  and  at  least  four  times  as  skillful  as 

the  best  operational  model  in  the  northwest  Pacific.  These  very 

encouraging  vertiflcation  results  suggest  the  potential  for  significant 

reduction  in  foreast  position  errors  and  should  hasten  the  testing  of 

the  model  in  an  operational  environment., 
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1.  INTRODUCTION 


The  problem  of  forecasting  the  movement  of  tropical  cyclones, 
particularly  recurvature  and  landfall,  is  crucial  in  many  parts  of  the 
world.  Considerable  effort  has  been  expended  by  the  various 
Meteorological  Services  of  a  number  of  nations  on  improving  the  accuracy 
of  track  prediction.  A  historical  overview  will  be  presented  in  Chapter 
2. 

Although  significant  progress  has  been  made  in  developing 
statistical  and  dynamical  models  for  forecasting  the  tracks  of  tropical 
cyclones,  the  forecast  errors  over  the  last  few  years  (Fig.  1.1)  have 
shown  little  improvement.  The  official  forecasts  for  the  NW  Pacific 
issued  by  the  US  Joint  Typhoon  Warning  Center  (JTWC)  on  Guam  have 
average  forecast  position  errors  on  the  order  of  120,  240,  and  360 
nautical  miles  (n  ml.)  for  24-,  48-,  and  72-hour  forecasts  while  those 
for  the  North  Atlantic  issued  by  the  US  National  Hurricane  Center  (NHC) 
in  Miami  have  errors  on  the  order  of  120,  250,  and  380  n  mi. 

Attempts  to  develop  more  sophisticated  numerical  models  oontinue  in 
hopes  of  achieving  a  breakthrough  in  track  prediction.  However,  the 
achievement  of  this  goal  does  not  appear  to  be  close  at  hand.  In 
addition,  in  most  developing  oountrles,  computer  facilities  are  not 
adequate  to  apply  the  more  sophisticated  dynamical  models.  Therefore, 
the  need  remains  for  a  tropical  cyclone  forecasting  technique  requiring 
only  meteorological  data  from  a  network  of  conventional  surface  and 
upper-air  observing  stations  as  basic  input. 
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ATLANTIC  FORECAST  ERRORS 
YEARLY  MEAN 


NW  PACIFIC  FORECAST  ERRORS 
YEARLY  MEAN 


Pig.  1.1.  The  yearly  seen  official  foreoaat  errors  (n  mi.)  during  the 
period  1974-1983  for  the  Atlantic  and  NW  Pacific.  The  yearly 
fluctuation  is  much  higher  in  the  Atlantic  than  the  NW 
Pacific  due  to  far  fewer  cases.  In  1977*  for  example,  there 
were  only  two  cases  for  verification  at  72  hours.  The 
1974-1982  errors  for  the  Atlantic  are  from  Annual  Data  and 
Verification  Tabulation  while  the  1983  errors  are  from  the 
minutes  of  the  38th  Annual  International  Hurricane  Conference 
The  errors  for  the  NW  Pacific  are  from  the  Annual  Tropical 
Cyclone  Typhoon  Reports.  1974-1983. 
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Furthermore,  a  review  of  objective  techniques  currently  in 
operational  use  reveals  inherent  weaknesses  in  their  design.  These 
weaknesses  will  be  specifically  addressed  in  Chapter  3. 

The  purpose  of  this  research  was  to  apply  a  new  statistical  method 
of  introducing  synoptic  and  persistence  predictors  to  develop  an 
Improved  prediction  model  that  could  be  used  at  all  operational  forecast 
centers.  Such  a  model  would  need  to  be  simple  so  that  it  could  operate 
without  a  link  to  a  large  computer  data  base  or  numerically  produced 
analysis  or  forecast  fields,  yet  be  sophisticated  enough  to  produce 
accurate  forecasts  at  both  the  short  and  longer  forecast  periods. 

The  motivation  for  this  work  arose  as  a  result  of  nearly  three 
years  of  operational  tropical  cyclone  forecasting  by  the  author  during 
which  time  the  strengths  and  weaknesses  of  the  myriad  of  objective 
techniques  in  use  became  obvious.  Almost  all  of  these  track  prediction 
schemes  perform  adequately  on  well-behaved  cyclones  that  follow 
climatological  tracks.  However,  on  those  anomalous,  difficult  cases 
when  the  forecaster  needs  objective  guidance  the  most,  the  schemes 
Invariably  fail. 

In  developing  a  new  and  improved  statistical  forecast  model,  the 
method  of  Introducing  and  applying  the  various  climatological, 
persistence  and  synoptic  predictors  for  predicting  the  future  motion  of 
tropical  cyclones  is  patterned  after  the  normal  procedure  followed  by 
most  forecasters  as  they  subjectively  formulate  the  forecast  track. 

The  innovative  approach  developed  here  primarily  Involves  the 
discretizing  of  the  72-hour  foreoast  track  into  three  24-hour  time  steps 
and  applying  updated  persistence  and  synoptic  data  relative  to  each  new 
cyclone  position  as  the  forecast  is  stepped  forward  in  24-hour 


increments.  The  perfect-prog  method  is  used  to  introduce  the  updated 
synoptic  Information  but  only  at  the  valid  time  of  the  prognostic  data 
and  always  relative  to  the  cyclone  position  at  the  valid  time  of  the 
prognostic  data.  Climatology  is  incorporated  via  a  stratification 
scheme  based  on  the  position  of  the  cyclone  relative  to  the  500  mb 
subtropical  ridge.  Segmenting  the  72-hour  cyclone  motion  permits  the 
category  the  cyclone  is  in  to  change  after  each  24  hours  of  motion,  if 
necessary.  The  three  categories  relative  to  the  ridge  —  south,  on, 
north  —  result  in  three  sets  of  regression  equations  with  one  set 
comprised  of  equations  for  the  north-south  and  east-west  displacements. 

A  more  detailed  discussion  of  the  model  designs  will  be  presented  in 
Chapter  3. 

Chapter  4  will  discuss  the  multiple  regression  technique  used  to 
handle  the  large  number  of  cases  and  potential  predictors.  Correlation 
coefficient  fields  between  the  synoptic  predictors  and  east-west  and 
north-south  cyclone  motion  will  be  presented  along  with  an  examination 
of  the  contribution  of  each  of  the  predictors  to  the  reduction  of 
variance. 

Results  from  the  verification  of  this  new  statistical  prediction 
model,  to  be  referred  to  as  CSU84,  will  be  presented  in  Chapter  5.  Also 
presented  in  this  chapter  will  be  sensitivity  tests  that  simulate 
operational  conditions.  This  was  accomplished  by  perturbing  the  "best 
track  initial  positions  and  by  using  actual  forecast  fields 
produced  by  a  numerical  model. 

***  The  best  track  positions  are  the  accepted  cyclone  positions 

after  a  post-cyclone  analysis. 


Although  CSU84  was  initially  an  attempt  to  improve  track  forecasts 
primarily  at  the  higher  latitudes  and  at  the  longer  forecast  periods, 
the  results  indicate  considerable  skill  at  the  lower  latitudes  and 
shorter  forecast  periods  as  well. 


2.  BACKGROUND 


The  movements  of  tropical  cyclones  remain  one  of  the  most  important 
and  difficult  forecasting  problems  in  tropical  meteorology.  Earliest 
attempts  at  developing  objective  forecast  schemes  hypothesized  that 
tropical  cyclones  were  carried  along  by  the  prevailing  larger  scale 
environmental  flows  in  which  they  were  embedded  like  eddies  in  a  river. 
This  concept  has  come  to  be  known  as  "steering". 

More  recently,  observational  studies  have  demonstrated  a  highly 
significant  relationship  between  tropical  cyclone  motion  and  synoptic- 
scale  flow  (**  500-1000  km)  around  the  cyclone.  The  pressure  level  at 
which  the  speed  and  direction  of  the  surrounding  winds  (or  equivalently 
the  pressure  or  height  gradients  across  the  cyclone)  best  correlate  with 
those  of  the  cyclone  has  been  referred  to  as  the  steering  level. 

One  problem  in  applying  the  steering  concept  has  been  the  selection 
of  the  proper  level  that  best  correlated  with  the  movement  of  the 
cyclone.  Although  different  forecast  schemes  use  different  steering 
levels,  it  is  generally  accepted  that  the  mid-tropospheric  levels  are 
the  best  in  predicting  tropical  cyclone  movement. 

George  (1975)  found  that  for  NW  Pacific  cyclones,  the  composite  500 
mb  winds  averaged  between  1-7°  latitude  radius  from  the  cyclone  center 
had  the  strongest  correlation  with  the  direction  of  cyclone  movement 
while  the  700  mb  winds  best  correlated  with  the  cyclone  speed.  Gray 
(1977)  presented  similar  results  for  Atlantic  oyclones.  Chan  (1982) 
using  composite  wind  data  over  an  area  (5-7°  latitude  radius  from  the 


cyclone  center)  outside  the  strong  inner  circulation  of  the  cyclone 
found  that  the  winds  at  the  aid-troposphere  (700,  600  and  500  mb) 
correlated  best  with  both  the  direction  and  speed  of  cyclone  movement. 
More  importantly,  Chan  found  that  cyclones  having  different  zonal 
components  of  motion  have  different  relationships  with  their  5-7° 
surrounding  flow.  Holland  (1983)  showed  that  cyclone  motion  was  the 
result  of  a  nonlinear  combination  of  two  processes:  1)  an  Interaction 
between  the  cyclone  and  its  basic  current  (the  well-known  steering 
concept),  and  2)  an  interaction  with  the  earth's  vorticity  field  which 
causes  a  westward  deviation  from  the  pure  steering  flow.  Similar  to 
Chan's  findings,  Holland  discovered  that  westward  moving  cyclones  were 
only  marginally  affected  by  fluctuations  in  the  basic  flow  whereas 
northward  and  eastward  moving  cyclones  may  undergo  large  motion  changes 
for  a  similar  range  of  fluctuations. 

Attempts  to  use  winds  and  heights  at  upper  tropospheric  levels 
(e.g..  Jordan,  1952;  Miller,  1958)  have  historically  not  been  as 
successful  and  have  received  little  use.  However,  recent  findings  have 
supported  a  good  correlation  between  upper  tropospheric  winds  and  the 
recurvature  phenomenon.  Specifically,  George  and  Gray  (1977),  using  10 
years  of  composited  rawinsonde  reports,  found  large-scale  differences  at 
the  200  mb  level  15°  to  20°  from  the  cyclone  prior  to  significant 
changes  in  cyclone  movement.  Their  results  indicated  that  the  strength 
of  the  upper  level  westerlies  northwest  and  north  of  the  cyclone  is 
strongly  related  to  whether  a  cyclone  recurved  or  not. 

Bao  and  Sadler  (1982)  also  found  a  strong  relationship  between  the 
200  mb  winds  at  and  before  recurvature  along  the  cyclone's  subsequent 
track  and  the  speed  of  movement  of  the  cyclone  after  recurvature. 
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Gentry  (1983)  found  that  when  the  mean  200  mb  winds  at  about  15°  to 
20°  northwest  and  north  of  the  cyclone  equal  or  exceed  20  m/s  that  80% 
of  the  cyclones  recurved  before  traveling  as  much  as  12  degrees  of 
longitude  further  west. 

Weir  (1982)  developed  an  operational  technique  for  the  NW  Pacific 
using  200  mb  winds  to  assist  in  forecasting  the  acceleration  of  typhoons 
as  they  Interact  with  the  mid-latitude  westerlies. 

Apparently,  upper  tropospheric  levels  are  generally  poor  steering 
levels  with  the  exception  of  recurving  cyclones.  The  reason  for  this 
observed  peculiarity  is  perhaps  due  to  the  greater  200-500  mb 
directional  shear  in  the  deep  tropics.  With  recurving  cyclones,  on  the 
other  hand,  directional  shear  is  much  less  and  the  200  mb  winds  provide 
an  extra  degree  of  predictive  Information  for  cyclones  interacting  with 
the  mid-latitude  westerlies. 

The  surface-pressure  map  together  with  24-hour  surface-pressure 
changes  (to  eliminate  large  diurnal  effects  in  the  tropics)  has  been 
used  by  tropical  cyclone  forecast  centers  since  the  1950 's.  It  is  still 
used  in  conjunction  with  more  sophisticated  techniques  to  better  assess 
the  shorter  range  objective  forecasts  versus  the  current  synoptic 
situation.  In  addition,  pressure  changes  outside  the  cyclone's 
circulation  have  been  used  to  indicate  future  changes  in  motion.  When 
there  are  large  pressure  changes  far  removed  from  the  cyclone,  their 
effect  is  to  turn  the  cyclone  at  right  angles  to  the  line  connecting  the 
isallobaric  center  and  the  cyclone  center. 

A  new  prediction  model  was  developed  with  the  purpose  of 
assimilating  all  of  these  observations  linking  cyclone  motion  to  various 
synoptic  scale  flow.  Specifically,  all  cyclones  were  stratified 


relative  to  the  500  mb  subtropical  ridge  to  incorporate  the  findings  of 
Chan  (1982),  Xu  and  Gray  (1982).  Holland  (1983)  and  others  indicating 
different  steering  relationships  for  cyclones  having  different  zonal 
components  of  motion  or  different  locations  relative  to  the  ridge. 
Further,  200  mb  steering  predictors  were  used  0-20°  northwest  and  north 
of  the  cyclone  based  on  the  correlation  between  upper  tropospheric  winds 
and  recurvature  shown  by  the  numerous  studies  just  cited.  Five  hundred 
millibar  heights  and  surface  pressures  and  their  24-hour  changes  were 
used  as  predictors  to  assess  the  current  and  recent  changes  in  the 
synoptic  situation  at  two  different  levels.  The  use  of  200  mb  steering 
winds  provided  a  third  vertical  level.  Seventy-two  hours  of  cyclone 
motion  were  segmented  into  three  shorter  24-hour  time  steps  and  forecast 
fields  relative  to  each  new  cyclone  position  were  used  to  better  define 
the  synoptic  flow  and  its  changes  affecting  the  movement  of  the  cyclone 
as  it  evolved  over  a  3-day  period. 

Based  on  the  steering  concept  just  mentioned,  a  number  of  other 
tropical  cyclone  track  forecasting  schemes  have  been  developed  to  relate 
statistically  or  dynamically  cyclone  motion  to  the  surrounding  flow 
field.  A  review  of  the  basic  categories  of  objective  techniques  used 
for  predicting  tropical  cyclone  motion  and  a  general  summary  of  the 
major  tropical  cyclone  track  prediction  methods  developed  over  the  years 
will  now  be  presented. 

There  are  four  basic  categories  of  objective  techniques  used  for 
predicting  tropical  cyclone  motion: 

(1)  climatology-persistence 

(2)  steering  flow 

(3)  statistical 

(4)  numerical 


Appendix  A  presents  the  various  techniques  currently  in  operational 
use  in  the  Atlantic  and  the  NW  Pacific/North  Indian  Oceans.  A  more 
detailed  discussion  of  these  track  prediction  models  can  be  found  in 
Neumann  and  Peilissier  (1981)  for  the  Atlantic  and  in  the  Annual 
Tropical  Cyclone  Report  prepared  by  the  Joint  Xyphoon  gaming  (Renter 
(JTWC)  for  the  other  ocean  basins. 

2.1  Climatology-Persistence  Methods 
Climatology-persistence  methods  rely  upon  empirical  relationships 

related  to  the  tracks  of  previous  cyclones  or  the  past  motion  of  the 
current  cyclone.  These  techniques  are  useful  for  well-behaved  cyclones 
that  move  along  a  climatological  track  with  few  changes  in  direction. 
Because  most  cyclones  in  the  deep  tropics  embedded  in  the  easterlies 
exhibit  this  type  of  motion,  climatology-persistence  methods  perform 
well  for  cyclones  at  low  latitudes.  But  these  techniques  break  down  in 
the  westerly  wind  belt  regions. 

2.2  Steering  Flow  Method 

The  steering  flow  concept  views  the  tropical  cyclone  as  point 
vortices  embedded  in  the  basic  environmental  flow  and  being  moved  by  the 
steering  current.  Observations  support  this  ooncept  of  the  cyclone 
responding  to  the  environmental  flow  in  which  it  is  embedded. 

Difficulty  of  applying  this  concept  operationally  occurs  when  the  data 
does  not  allow  proper  steering  definition  or  when  vertical  shear  in  the 
environmental  flow  is  present  and  there  is  doubt  as  to  which  level  is 
the  optimum  steering  level. 

2.3  Statistical  Method 

The  statistical  approach  oompares  different  parameters  or 
combinations  of  parameters  quantitatively  to  cyclone  motion.  Multiple 


regression  analysis  is  employed  to  develop  correlation  coefficients  and 
regression  equations. 

As  summarized  by  Neumann  (1979),  statistical  models  for  predicting 
tropical  cyclone  motion  cam  be  grouped  into  two  broad  categories:  those 
models  based  on  analog  and  those  based  on  regression  equations.  The 
latter,  in  order  of  increasing  sophistication,  can  be  further 
categorized  as  follows: 

(a)  Those  models  which  use  predictors  based  on  climatology  and 
persistence: 

(b)  Those  models  which  include,  but  which  are  not  limited  to, 
predictors  derived  from  observed  synoptic  data: 

(c)  Those  models  which  include,  but  which  are  not  limited  to 
predictors  (a)  and  (b),  but  also  include  predictors  derived 
from  numerically-forecast  data. 

The  models  that  fall  into  category  (c)  are  known  as 
statistical/dynamical  models  while  those  that  fall  into  categories  (a) 
and  (b)  have  been  referred  to  respectively  as  simulated  analog  models 
and  classical  or  statistical-synoptic  models. 

2.3.1  Analog  Models 

2. 3. 1.1  The  Analog  Concept 

Analog  models  derive  their  success  from  their  ability  to  identify 
families  of  cyclone  tracks.  Using  a  series  of  computer  algorithms,  a 
current  cyclone  is  associated  with  its  most  likely  family,  thus  allowing 
inferences  to  be  made  about  its  future  motion.  The  analog  process 
identifies  historical  cyclones  temporally  and  spatially  similar  to  the 
current  cyclone.  Typically,  such  factors  as  time  of  year,  the  cyclone's 
initial  direction  of  motion  and  speed  and  the  position  of  the  analog 
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cyclone  relative  to  the  current  cyclone  are  considered.  Selected  analog 
tracks,  translated  to  a  common  origin,  and  combined  with  persistence, 
make  up  the  final  forecast.  Refinements  over  the  years  have  included 
analog  weighting  schemes,  use  of  past  as  well  as  current  motion  and 
better  combinations  of  climatology  and  persistence. 

2. 3. 1.2  Historical  Development 

The  first  fully  operational  model  for  predicting  tropical  cyclone 
motion  by  analog  methods  was  developed  by  Hope  and  Neumann  (1970)  for 
use  by  the  US  National  fiurricane  Center  (NHC)  on  Atlantic  cyclones.  The 
technique,  known  by  the  acronym  HURRAN  (HURRicane  ANalog).  has  been  in 
continuous  use  at  the  Center  since  the  beginning  of  the  1969  hurricane 
season. 

Concurrently  and  independently,  Hodge  and  McKay  (1970)  collaborated 
with  the  US  Navy  and  developed  an  analog  prediction  model  for  North 
Pacific  typhoons.  Their  model,  known  as  TYF00N,  was  revised  by  Jarrell 
and  Somerville  (1970)  and  became  operational  at  the  Joint  Typhoon 
Warning  Center,  Guam  in  August  1970. 

Following  the  introduction  of  HURRAN  and  TYF00N,  and  their  moderate 
success,  the  analog  concept  became  quite  popular.  The  HURRAN/TYF00N 
system  was  adopted  and  modified  for  use  in  the  North  Indian  Ocean, 
southwest  Indian  Ocean,  northeastern  Pacific  and  in  the  southwest 
Pacific  and  Australian  area. 

2. 3. 1.3  Advantages  and  Disadvantages  of  the  Analog  Method 

The  main  advantage  of  the  analog  is  its  simplicity  and  economy  in 
terms  of  computer  resources.  It  is  quick  to  run  and  is  one  of  the  first 
objective  guidances  available  to  the  forecaster.  In  addition,  the 
method  is  very  accurate  for  cyclones  in  the  deep  tropics  embedded  in  the 


easterlies  since  analog  forecasts  represent  typical  or  normal  cyclone 
motion  and  deviations  from  normal  in  the  tropics  are  small. 

However,  at  higher  latitudes  where  cyclones  are  affected  by  the 
westerlies,  anomalous  situations  are  much  more  oommon  and  analog 
forecasts  are  more  unreliable.  Neumann  and  Hope  (1972)  have  found  that 
about  one  out  of  three  forecast  situations  in  the  Atlantic  is  too 
anomalous  for  an  analog  method  to  arrive  at  a  solution. 

Another  disadvantage  of  the  analog  method  is  a  forecast  bias 
Introduced  when  cyclones  are  prematurely  dropped  from  the  history  tape 
due  to  landfall,  extratropieal  transition  or  dissipation.  In  these 
instances,  the  forecast  is  excessively  influenced  by  the  remaining 
cyclone  tracks. 

2.3.2  Regression  Equation  Models 
2. 3. 2.1  Historical  Development 

The  first  of  the  statistical  systems  for  tropical  cyclone  track 
forecasting  was  developed  by  Riehl,  .§£  al.  (1956)  using  the  geostrophic 
wind  components  computed  from  the  500  mb  heights  around  the  periphery  of 
the  cyclone  as  a  method  of  representing  the  vertically  integrated 
steering  flow.  Miller  and  Moore  (1960)  calculated  additional 
geostrophic  wind  components  for  steering  at  the  700  and  300  mb  levels  as 
well  as  the  past  movement  of  the  cyclone  as  predictors.  Arakawa  (1964) 
and  Tse  (1966)  developed  a  similar  method  for  use  in  the  Pacific  as  did 
Kumar  and  Prasad  (1973)  for  the  North  Indian  Ocean.  The  Tse  method  is 
somewhat  unique  in  that  one  of  five  sets  of  regression  equations  were 
used,  depending  on  which  one  of  five  easily  Identifiable  700  mb  patterns 


existed. 


All  of  these  methods  used  a  relatively  small  grid  surrounding  the 
cyclone  to  forecast  24-hour  cyclone  motion.  However,  because  of  their 
sensitivity  to  the  initial  analysis,  all  suffer  from  the  subjectivity  of 
uncertain  analysis  close  to  the  cyclone  center. 

Veigas,  £&  (1959)  were  the  first  to  apply  screening  -  multiple 

linear  regression  statistical  methods  to  the  problem  of  tropical  cyclone 
forecasting.  Rather  than  restrict  potential  predictors  to  the  cyclone 
vicinity.  Veigas  used  a  large  storm-centered  grid  to  introduce 
predictors  that  represented  large-scale  weather  patterns.  This  method 
had  the  advantage  of  being  less  sensitive  to  the  initial  analysis  near 
the  cyclone. 

A  major  advance  in  statistical  prediction  of  tropical  cyclone 
motion  occurred  with  the  development  of  the  NHC64  system.  In  this 
scheme.  Miller  and  Chase  (1966)  combined  the  better  features  of  the 
Veigas.  Mlller-Moore  and  Rlehl-Haggard  methods  into  a  single  model. 

Sea  level  pressures.  700  mb  heights  and  500  mb  heights  as  well  as 
1000-700  mb  thicknesses.  700-500  mb  thicknesses,  500  mb  height  changes, 
geostrophic  wind  components  at  the  1000  mb,  700  mb  and  500  mb  levels  and 
the  past  12-hour  movement  of  the  cyclone  center  were  used  as  predictors. 
Cyclones  greater  than  34  kts  for  the  period  1945-1961  formed  the 
development  set  which  was  divided  into  two  classes  -  north  and  south  of 
27.5°N.  A  screening-multiple  linear  regression  technique  was  used  to 
derive  four  sets  of  equations  —  one  set  for  the  0-12  hour  north-south 
and  east-west  displacements,  a  second  set  for  the  12-24  hour 

displacements,  a  third  set  for  the  24-36  hour  displacements  and  a  fourth 

% 

set  for  the  36-48  hour  displacements.  Forecasts  for  periods  greater 
than  12  hours  were  prepared  by  adding  forecast  displacements  for 


successive  periods.  A  48-72  hour  forecast  was  added  to  the  technique 
later. 

As  a  result  of  NHC64's  poor  performance  on  an  erratic  cyclone 
(Betsy)  in  1966,  a  new  statistical  model  was  developed  (Miller,  si  al. . 
1968).  The  revised  model,  known  as  NHC67,  is  still  in  use  at  the 
National  Hurricane  Center  in  Miami.  Data  from  hurricanes  during  the 
years  1962-1965  were  added  to  the  1945-1961  period  used  in  deriving  the 
NHC64  equations.  Additional  predictors  were  screened  in  developing  the 
NHC67  equations.  These  additional  predictors  were  the  1000  mb  heights 
(in  place  of  sea  level  pressure),  1000  and  700  mb  24-hour  height  changes 
and  the  thermal  winds  computed  from  the  thicknesses.  Although  equations 
to  make  the  48-72  hour  forecast  were  derived  due  to  the  operational 
requirement  to  issue  72-hour  forecasts,  it  was  felt  that  these  equations 
were  of  questionable  value. 

In  order  to  provide  an  analog  forecast  as  input  to  a  higher  echelon 
prediction  model  quickly  and  under  any  anomalous  synoptic  condition, 
Neumann  (1972)  developed  a  rather  simple  regression  equation  model  for 
the  Atlantic  area  known  by  the  acronym  CLIPER  (CLImatology  and 
PERslstence) .  Operational  use  of  this  improved  analog  model  confirmed 
that  climatology  and  persistence  were  being  better  used  than  by  earlier 
models  for  cyclones  embedded  in  the  easterlies.  However,  when  cyclones 
recurved  out  of  the  tropics,  the  analog  models  were  still  inferior  to 
older  models  such  as  NHC67  which  were  cognizant  of  environmental  flow 
patterns.  Accordingly,  Neumann,  si-  (1972)  blended  the  better 
features  of  both  approaches  into  a  new  model  for  the  Atlantic  area  known 


as  NHC72. 
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However,  Atlantic  cyclones  that  exhibited  anomalous  motion 
continued  to  confirm  the  inherent  inability  of  the  purely  "classical*' 
models  typified  by  NHC67  and  NHC72  to  forecast  such  motion.  This 
finally  resulted  in  the  development  of  the  so-called 
statist! cal/dynamical  (as  well  as  the  purely  dynamical)  models  which 
conceptually  should  respond  better  to  anomalous  synoptic  conditions. 

Such  models  use  the  output  from  a  numerical  model  as  one  set  of 
predictors  for  input.  Neumann  and  Lawrence  (1973)  used  selected 
predictors  from  the  24-,  36-,  and  48-hour  prognostic  500  mb  geopotential 
height  fields  produced  by  the  US  National  Meteorological  Center  (NMC) 
Primitive  Equation  (PE)  model  for  their  statistical/dynamical  NHC73 
model.  A  great  deal  of  success  was  achieved  by  NHC73,  primarily  on  the 
longer-range  forecasts.  Another  statistical/dynamical  model  developed 
by  Nomoto,  ifc  si-  (1976)  is  being  used  by  the  Japan  Meteorological 
Agency  on  NW  Pacific  typhoons  in  the  vicinity  of  Japan. 

The  term  statistical/dynamical  can  also  be  applied  to  the  HATTRACK 
(Hurricane  And  Typhoon- TRACKing)  (Renard,  1968)  and  the  MOHATT  (Modified 
HATTrack)  (Renard.  jgfc  al. ■  1973)  techniques  used  by  JTW C  on  Guam. 

Instead  of  regression  equations,  tropical  cyclones  are  geostrophlcally 
steered  in  six-hourly  time  steps  using  the  US  Navy's  smoothed  numerical 
prognoses  at  various  levels  (1000  mb,  850  mb,  700  mb,  500  mb,  400  mb, 
and  200  mb).  The  operational  forecasts  starting  at  time  Tq  are  based  on 
a  forecast  cycle  starting  at  time  Tq  minus  12  hours.  This  allows  for  an 
error  measurement  of  the  motion  for  the  first  12  hours  of  an  84-hour 
projection.  This  error,  assumed  to  be  a  bias,  is  linearly  applied  with 
certain  constraints,  to  the  remaining  72-bour  forecast  track.  The 


operational  version  of  HATTBACK/MOHATT  is  known  by  the  acronym  CYCLOPS 
(CYCLone  Operational  prediction  System) . 

Since  1983,  a  model  output  statistic  technique  known  by  the  acronym 
COSMOS  (£yclops  Objective  Steering  Model  Output  Statistics)  (Allen, 

1984)  has  been  in  operational  use  at  JTWC,  Guam,  to  assist  the 
forecaster  in  Interpreting  the  output  from  CYCLOPS.  COSMOS  uses  the 
CYCLOPS  forecasts  at  the  850,  700,  and  5 00  mb  levels  to  produce  its  own 
forecast  based  on  a  statistical  analysis  of  the  past  performance  of 
CYCLOPS. 

2. 3. 2. 2  Advantages  and  Disadvantages  of  the  "Classical"  Models 

Several  models  in  this  category  are  in  current  operational  use  at 
the  various  tropical  cyclone  forecast  oenters.  These  include  Arakawa’s 
(1964)  and  Tse's  (1966)  models  in  the  NW  Pacific;  Kumar  and  Prasad's 
(1973)  model  in  the  North  Indian  Ocean;  and  the  North  Atlantic  NHC67 
(Miller,  fiifil. ,  1968)  and  NHC72  (Neumann,  fit  *1. ,  1972)  models.  All  of 
these  models  incorporate  predictors  obtained  from  current  or  recently 
observed  synoptic  data  to  forecast  future  positions  of  the  cyclone  out 
to  as  much  as  72  hours. 

The  reduction  of  variance  of  such  models  typically  decreases 
rapidly  with  increased  forecast  projection.  Predictors  are  selected 
from  climatology  and  persistence  as  well  as  one  or  more  of  the  current 
and  24-hour  old  1000,  700,  and  500  mb  height  fields.  In  fact, 
persistence  provides  most  of  the  reduction  of  variance  for  the  shorter- 
term  forecast. 

Separate  sets  of  regression  equations  (stratification)  representing 
different  synoptic  patterns  have  been  developed  for  some  of  these 
models.  Tse's  (1966)  model  stratifies  according  to  five  easily 


Identifiable  700  mb  patterns.  The  NHC67  nod el  stratifies  according  to 
the  cyclone's  initial  location  north  or  south  of  27.S°N  with  an 
additional  substratification  dependent  on  the  cyclone's  speed.  A 
stratification  based  on  the  cyclone's  initial  motion  is  used  by  the 
NHC72  model. 

Both  the  NHC67  and  the  NHC72  (as  well  as  the  NHC73)  models  use 

predictors  from  three  levels  —  1000,  700,  and  500  mb.  According  to 

Neumann  (1979),  tests  conducted  at  the  NHC  subsequent  to  the  development 

of  these  models  suggest  that  better  operational  performance  may  be 

possible  by  the  elimination  of  the  700  mb  data.  Data  at  this  level  are 

highly  correlated  with  the  500  mb  level  and  it  appears  that  the  models 

use  more  predictors  than  warranted  by  the  sample  size. 

2. 3. 2. 3  Advantages  and  Disadvantages  of  the  Statistical- 
Dynamical  Models 

At  the  top  of  the  statistical  models  for  the  prediction  of  tropical 
cyclone  motion  are  those  models  which  derive  a  portion  of  the  variance 
reduction  from  predictors  taken  from  the  output  of  a  numerical 
prediction  model.  Such  models  also  derive  a  portion  of  their  forecast 
skill  from  other  sources  such  as  climatology,  persistence  or  current 
synoptic  data. 

The  use  of  jjodel  jQutput  Statistics  (M0S)  is  conceptually  the  most 
effective  method  of  introducing  numerically  forecast  data  into  a 
statistical  prediction  model.  However,  the  lack  of  adequate  data  sample 
prevents  the  direct  Introduction  of  a  numerical  model  output  into  a 
statistical  model.  As  a  result,  substitute  methods  of  using  the 
statistical -dynamical  approach  have  been  derived.  One  of  these  is 
referred  to  as  the  "perfect-prog"  approach  and  a  second  as  the 


Simulated  Model  fiutput  Statistics  (SMOS).  The  perfect-prog  approach  was 
used  in  the  development  of  the  model  described  in  this  paper. 

In  the  perfect-prog  method,  observed  values  of  a  predictor  at  time 
TQ  +  AT  are  used  to  derive  a  statistical  relationship  between  predictor 
and  predictand  at  the  same  time  Tq  +  AT.  In  actual  practice,  forecast 
values  of  the  predictor  must  be  used  and,  therefore,  any  forecast  error 
or  bias  is  passed  on  to  the  statistical  model.  This  is  one  reason  why 
there  is  a  tendency  for  the  predictive  skill  of  statist! cal -dynamical 
models  developed  using  the  perfect-prog  method  to  degrade  significantly 
when  used  operationally.  With  the  MOS  technique,  the  biases  of  the 
numerical  model  would  have  been  statistically  corrected  and  certain 
inaccuracies  recognized. 

In  the  SHOS  method,  the  covariance  matrices  containing  the  sums  and 
cross-products  of  all  possible  combinations  of  variable  pairs  are 
modified  or  "contaminated"  by  introducing  an  error  component  similar 
to  that  of  actual  prognostic  data.  The  purpose  of  the  contamination  is 
to  partially  correct  the  tendency  for  the  screening  program  to  over¬ 
weight  or  select  too  many  predictors  from  the  perfect-prog  data  at  the 
expense  of  the  other  potential  predictors. 

As  discussed  in  Neumann  and  Lawrence  (1973),  NHC73  was  tested  on 
operational  independent  data  using  different  sets  of  regression 
equations  developed  from  the  two  methods.  Both  versions  performed  well 
at  the  extended  forecast  periods.  As  expected,  the  perfect-prog 
equations  performed  relatively  poorly  during  the  first  half  of  the  72- 
hour  forecast  period.  Predictor  over-weighting  was  assumed  to  be  the 
probable  explanation.  The  Improvement  in  the  perfect-prog  forecasts  at 
the  extended  periods  suggested  that  the  variance  reducing  potential  of 
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the  numerically  forecast  data  Is  greater  than  that  provided  by  the 
current  data  at  48  hours  and  beyond. 

The  better  operational  performance  of  the  SMOS  equations  for  the 
shorter  range  forecast  periods  was  believed  to  be  the  more  realistic 
blend  of  current  and  numerically  forecast  data  as  a  result  of  the 
contamination  of  the  perfect-prog  predictors. 

In  Japan,  a  statistical -dynamical  model  described  by  Nomoto,  £&  al. 
(1976)  is  also  designed  around  the  perfect-prog  concept.  Its  synoptic 
predictors  oome  only  ft*om  numerically  produced  24-hour  500  mb  prognositc 
geopotential  height  fields.  With  additional  persistence  predictors,  the 
model  provides  12-.  24-.  36-.  and  48-hour  forecast  positions.  Vith  the 
exception  of  the  HATTRACK/MOHATT  version  called  CYCLOPS,  there  are  no 


Numerical  models  "appear"  to  provide  the  best  promise  for 
improvement  in  track  forecasting.  Dynamical  techniques  assume  that  the 
motion  of  the  cyclone  may  be  forecast  directly  from  numerical 
Integration  of  geophysical  governing  equations  such  as  momentum, 
continuity  and  thermodynamics. 

2.4.1  Historical  Development 

Harrison  (1973)  developed  a  simple  nested  grid  model  to  forecast 
typhoon  motion  in  the  NW  Pacific  using  the  primitive  equations.  This  is 
the  original  version  of  the  nested  two-way  Interactive  tropical  cyclone 
model  (NTCM)  in  operational  use  at  JTWC  on  Guam  (Harrison,  1981).  Prior 
to  the  development  of  the  operational  version  of  NTCM,  a  more  simple 
version  of  the  model  without  the  fine  grid  was  adapted  for  use  in  1975. 


This  coarse-mesh  (2°).  three-layer  dry  model  initially  possessed 
boundary  conditions  which  were  insulated,  free  slip  walls  on  the  north 
and  south  and  cyclic  on  the  east  and  west.  Hodur  and  Burk  (1978) 
modified  the  "channel"  model  by  incorporating  one-way  Interactive 
boundaries  with  forecast  fields  from  a  global/hemispheric  PE  model 
providing  the  boundary  values  in  12  hour  intervals.  This  modification, 
along  with  the  adjustment  of  the  initial  wind  fields  near  the  center 
based  on  the  difference  between  the  actual  and  predicted  cyclone  motion 
after  6  hours  using  a  technique  developed  by  Shewchuk  and  Elsberry 
(1978).  are  Incorporated  in  the  current  operational  model  known  as  OTCM. 
The  NTCM  and  the  OTCM  are  two  of  the  better  objective  techniques  at  JTWC 
providing  fairly  accurate  overall  track  prediction. 

In  the  Atlantic,  the  SANBAR  barotropic  model  (Sanders  and  Burpee, 
1968)  was  developed  on  the  belief  that  momentum  advection  is  the  primary 
physical  mechanism  for  motion  of  Intense  tropical  cyclones  and  that  the 
cyclone  is  steered  by  the  large-scale  environmental  flow  in  which  it  is 
embedded.  Accordingly,  the  model  uses  a  deep-layer  (1000  -  100  mb) 
pressure  -  weighted  U  and  V  wind  field  analysis  over  the  grid  domain  of 
the  model.  In  the  version  of  the  model  used  operationally  at  the  NHC,  a 
technique  described  by  Pike  (1972)  is  used  to  modify  the  wind  field  near 
the  cyclone  to  make  it  conform  better  to  the  initial  cyclone  motion 
vector.  This  technique  is  similar  to  one  developed  later  by  Shewchuk 
and  Elsberry  (1978)  for  use  on  the  OTCM  for  the  NW  Pacific. 

A  baroclinic  Movable  £ine  Mesh  (MFM)  dynamical  model  (Hovermale  and 
Livezey,  1977)  at  the  NHC  was  first  tested  operationally  on  Atlantic 
tropical  cyclones  during  the  1975  season  and  has  been  in  operational  use 
since  1976.  The  MFM  has  far  more  resolution  in  the  vertical  (10  layers) 


than  the  NTCM  or  the  OTCM  and  also  includes  parameterization  of  cumulus 
convection  and  the  planetary  boundary  layer  as  veil  as  topography.  A 
fine  grid  of  60  km  is  used  over  its  entire  domain. 

The  NTCM  and  OTCM.  on  the  other  hand,  have  only  three  vertical 
layers  and  use  a  simple  analytic  heating  scheme  with  no  topography.  The 
finer  scale  of  the  NTCM  (41  km  grid)  moves  with  the  cyclone  within  a 
coarser  channel  model  with  grid  spacing  of  205  km.  Once  initialized, 
the  model  runs  independently  of  any  forecast  fields.  By  contrast,  the 
boundary  conditions  of  both  the  MFM  and  OTCM  are  externally  specified  by 
forecast  fields  from  the  NMC  spectral  global  model  and  FNOC  spectral 
global  model,  respectively. 

2.4.2  Advantages  and  Disadvantages  of  Numerical  Models 

Since  a  dynamical  model  does  not  depend  on  the  statistical 
relationships  between  a  current  cyclone  and  historical  cyclones,  the 
reliability  of  the  model  could  be  expected  to  be  Independent  of  the 
details  of  a  cyclone's  past  behavior.  Thus,  variations  in  the  future 
track  of  a  cyclone  ought  to  be  Implicitly  predicted  through  the 
equations  of  motion. 

However,  due  to  uncertain  Initial  analyses,  principally  in  and 
around  the  cyclone  vortex,  all  dynamical  models  are  noted  for  poor 
performance  at  the  shorter  forecast  period  and  an  overall  motion  that  is 
too  slow.  They  do  demonstrate  increased  skill  for  the  longer  forecast 
periods  beyond  36  hours.  This  observation  has  been  noted  by  Florino,  si  &1 
(1982)  for  the  NTCM  and  Neumann  and  Pellssier  (1981)  for  the  MFM.  In  an 
effort  to  improve  the  accuracy  at  the  shorter  forecast  period,  both  the 
OTCM  and  SANBAft  incorporate  a  biasing  technique  that  forces  the  initial 
cyclone  motion  in  the  model  to  match  the  observed  motion  of  the  actual 


cyclone.  A  similar  bias-corrector  technique  is  envisioned  for  the  NTCM. 
Thus,  the  model  is  no  longer  independent  of  the  cyclone's  past  behavior 
which  is  touted  as  one  of  the  numerical  model's  strong  points. 

In  addition,  the  NTCM  is  initialized  from  the  operational  FNOC 
global  band  tropical  analysis  fields.  Since  the  analysis  contains  a 
mixture  of  observations,  observed  climatology  and  persistence,  in  data- 
sparce  areas  the  cyclone  track  will  probably  have  a  significant 
climatological  component. 

The  disadvantage  of  the  complex  and  resource  oostly  MFM  is  that 
computer  resource  limitations  preclude  any  more  than  a  48  hour  forecast, 
and  then  only  under  certain  high  threat  conditions.  Further,  guidance 
from  an  MFM  run  is  not  available  until  approximately  eight  hours  after 
synoptic  time.  Both  the  NTCM  and  OTCM,  due  to  their  relative 
simplicity,  are  run  routinely  to  72  hours  twice  a  day  on  all  tropical 
storms  and  typhoons.  Their  quick  turnaround  puts  the  model  output  at 
the  disposal  of  the  forecaster  at  the  same  time  as  he  receives  guidance 
from  the  other  objective  aids. 

2.5  The  Forecasting  Dilemma 

The  error  statistics  over  the  years  suggest  that  no  one  objective 
technique  performs  better  than  all  others  on  every  cyclone.  Some 
schemes  work  well  for  some  cyclones  but  poorly  on  others.  This  is  one 
of  the  reasons  for  the  myriad  of  forecast  aids  in  operational  use.  An 
obvious  problem  of  having  too  many  techniques  for  objective  guidance  is 
illustrated  in  Fig.  2.1  where  a  typical  "shotgun"  effect  makes  it 
difficult  for  the  forecaster  to  select  the  best  forecast  track.  Such  a 
scatter  in  the  forecast  track  is  common  when  a  large  directional  change, 
as  in  recurvature,  occurs.  Note  that  the  two  numerical  models  (NTCM  and 


OTCM)  and  COSMOS,  which  applies  MOS  technique  to  the  output  from  a 
purely  steering  scheme  at  three  vertical  levels,  performed  superior  to 
the  other  aids.  The  superior  performance  of  the  dynamical  models  over 
the  past  five  years  for  the  NW  Pacific  and  NHC73  for  the  Atlantic  north 
zone  indicates  the  importance  of  synoptic  data,  particularly  in  a 
prognostic  sense,  in  any  attempt  to  forecast  cyclone  motion. 

However,  the  lack  of  observational  data,  initialization  schemes 
that  disregard  certain  observational  data  in  favor  of  climatology,  and 
the  imperfect  physics  Involved  with  PE  models  combine  to  bring  about 
inaccurate  analysis  and  prognostic  fields.  If  the  forecaster  was  able 
to  modify  these  fields  when  they  were  obviously  in  error  an  Improved 


Fig.  2.1.  The  standard  array  of  JTWC’s  objective  forecasting  techniques 
available  to  support  the  19  August  1983  1200Z  warning  for 
Tropical  Storm  Dorn.  Included  are  the  forecast  issued  at 
191200Z  and  the  eventual  best  track.  (From  1983  Annual 
Tropical  Cyclone  Report.) 
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accuracy  in  track  prediction  would  likely  be  obtainable.  A 
statistical-dynamical  model  such  as  the  one  developed  and  detailed  in 
this  paper  permits  such  a  man-machine  interface. 

In  addition,  the  method  with  which  the  prognostic  data,  whether 
they  be  perfect-prog,  MOS  or  SMOS,  are  introduced  into  the  statistical 
scheme  appears  to  be  the  key  to  limiting  the  amount  of  degradation  all 
statistical-dynamical  models  suffer  when  they  are  used  operationally 
with  actual  forecast  fields.  An  Innovative  approach  to  incorporating 
the  prognostic  data,  as  well  as  persistence  and  climatology,  for  use  in 
a  statistical  model  has  been  developed  here  and  the  results  Indicate 
that  the  problem  of  operational  degradation  may  have  been,  to  a  large 
extent,  solved. 


3.  DEVELOPMENTAL  DATA  BASE  AND  METHODOLOGY 

Multiple  linear  regression  techniques  were  used  to  develop  forecast 
equations  for  the  prediction  of  the  east-west  (U)  and  north-south  (V) 
components  of  tropical  cyclone  motion.  Predictors  used,  as  shown  in 
Fig.  3.1,  were  past  12-  and  24-hour  motion,  steering  at  the  500  mb  and 
200  mb  levels  and  the  current  and  24-hour  changes  of  500  mb  geopotential 
heights  and  sea  level  pressures  relative  to  the  cyclone  position  at  each 
of  the  three  forecast  times  of  T=0,  T=24  hours,  T=48  hours.  This  key 
point  will  be  amplified  later. 

The  developmental  data  were  obtained  from  all  tropical  cyclone 
(including  tropical  depressions)  cases  between  1946-1981.  No  cases  were 
omitted  because  of  intensity  or  geographical  limitations.  The  year  1948 
and  every  third  year  thereafter  were  selected  to  be  used  as  an 
independent  set  for  testing.  Years  instead  of  cases  were  selected  for 
the  independent  set  because  it  was  felt  that  this  would  result  in  a  more 
representative  selection  for  all  three  categories  relative  to  the  ridge. 
Enough  years  were  selected  for  the  independent  set  to  provide 
approximately  one-third  the  number  of  cases  as  the  dependent  set.  The 
best  track  positions  for  the  Atlantic  basin  were  obtained  from  the  Best 
Track  of  The  National  Hurricane  Center  while  those  for  the  other  basins 
were  obtained  from  the  Annual  Typhoon /Tropical  Cyclone  Reports  published 
by  the  Joint  Typhoon  Warning  Center  (JTWC)  on  Guam.  The  synoptic  data 
were  derived  from  grid  point  analyses  archived  at  The  National  Center 
for  Atmospheric  Research  (NCAR).  These  grid  point  data  were 


Fig.  3.1.  The  three  components  of  predictors  used  by  CSU84. 


interpolated  onto  a  5  x  5  grid  (Fig.  3.2)  with  a  domain  of  -10  to 
+35°  latitude  and  ±40°  longitude  relative  to  the  cyclone  center  to 
represent  the  synoptic  forcing.  The  grid  was  offset  poleward  of  the 
cyclone  center  to  limit  the  number  of  grid  points  in  the  tropics,  where 
the  analyses  tend  to  revert  to  climatology.  A  similar  grid  has  been 
previously  used  by  Xu  and  Gray  (1982). 
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Fig.  3.2.  S°xS°  geographical  grid  with  a  domain  of  -10  to  +  35°  latitude 
and  +40°  longitude  relative  to  the  cyclone  center. 

Actual  prognostic  fields  of  S00  mb  D  values  and  sea  level  pressures 
used  to  verily  the  model  in  a  simulated  operational  setting  were 
received  from  the  Fleet  Numerical  Oceanography  Center  (FNOC),  Monterey, 
CA.  archived  on  a  63  x  63  hemispheric  grid  for  the  period  1968-1981. 

The  200  mb  prognostic  heights  were  not  available  from  the  Navy.  Thus, 
testing  of  the  model  using  actual  progs  was  limited  to  a  version  of  the 
model  that  did  not  include  200  mb  steering  predictors.  The  lack  of  200 
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mb  steering  predictors  did  not  present  a  problem  for  the  NW  Pacific 
since  testing  of  the  two  versions  revealed  no  significant  differences  in 
that  ocean  basin.  However,  differences  in  the  Atlantic  were  sizeable 
and  will  be  discussed  in  Chapter  S. 

3.1  Stratification  Scheme 

Thirty-six  years  of  tropical  cyclone  tracks  (1946-1981)  were 
subjectively  analyzed  along  with  daily  500  mb  synoptic  charts  to  select 
the  recurvature  point  for  each  cyclone.  Based  on  the  recurvature  point, 
every  00Z  and  12Z  position  along  with  its  subsequent  24-hour  track 
segment  was  placed  in  one  of  three  categories  relative  to  the  500  mb 
subtropical  ridge. 

The  south-of-ridge  category  is  oomprised  of  all  cyclones  that  never 
recurved  and  the  track  segment  of  those  recurving  cyclones  up  to  24 
hours  prior  to  the  recurvature  point.  The  north-of-ridge  category  is 
comprised  of  all  cyclones  that  existed  north  of  the  ridge  throughout 
their  lifetime  and  the  track  segment  of  those  recurving  cyclones  beyond 
24  hours  after  the  recurvature  point.  The  on-the-ridge  category  is 
comprised  of  the  track  segment  of  recurving  cyclones  between  24  hours 
before  and  24  hours  after  the  recurvature  point. 

Because  of  the  frequent  observation  of  recurving  cyclones  in  the  NW 
Pacific  that  exhibit  long  northward  tracks  before  turning  eastward,  the 
criteria  for  stratifying  the  cyclone  tracks  relative  to  the  ridge  were 
modified  for  the  NW  Pacific  and  North  Indian  Oceans.  Instead,  the 
south-of-ridge  category  includes  those  recurving  cyclones  up  to  the 
point  where  their  direction  of  motion  (DOM)  became  greater  than  330°. 

The  north-of-ridge  category  includes  those  recurving  cyclones  beyond  the 
point  where  their  DOM  became  greater  than  30°.  The  on-the-ridge 


category  is  comprised  of  the  track  segment  of  recurving  cyclones  where 
the  DON  is  between  330°  and  30°. 

In  an  attempt  to  better  forecast  the  extreme  cases  (which 
statistical  methods  frequently  fail  to  do),  the  three  categories  were 
further  separated  into  fast  and  slow  motion  groups  based  on  the  long 
wave  pattern  as  depicted  by  the  500  mb  height  difference  at  certain 
longitudes  along  a  chosen  latitude.  The  specific  grid  points  at  which 
the  height  differences  were  noted  were  determined  from  the  composition 
of  fast  minus  slow  cyclones  performed  by  Xu  and  Gray  (1982)  as  shown  in 
Figs.  3.3  and  3.4  for  the  Atlantic  and  NW  Pacific.  Based  on  their 
composition,  a  quantity  to  be  used  for  categorizing  the  speed  for  the 
Atlantic  cyclones  was  calculated  according  to: 

h15,-40  +  h15.25  “  2  X  Hij#-5  for  cyclones  north  of  the  ridge 

H10,10  "  h10.-15  "  20  for  cyclones  on  the  ridge 

2  *  h20,5  ~  h20,-25  ~  h20,25  for  cyclones  south  of  the  ridge 

and  for  the  Pacific: 

H10,-35  +  H10  is  “  2  x  Hio  -5  for  oydones  north  of  the  ridge 

H10  10  “  H10  -5  for  cyclones  on  the  ridge 

2  x  H15,0  “  H15,-30  ”  H15,40  for  cyclones  south  of  the  ridge 

H  is  the  S00  mb  height  in  meters  and  the  subscripts  indicate  the 
degrees  of  latitude  and  longitude  away  from  the  cyclone  center  where  the 
heights  are  to  be  observed.  Latitudes  are  positive  poleward  of  the 
cyclone  and  longitudes  are  positive  eastward  of  the  cyclone  center. 

If  the  calculated  quantity  was  greater  than  or  equal  to  zero,  the 
long  wave  pattern  was  assumed  to  be  favorable  for  fast  motion  so  the 
cyclone  was  placed  in  the  fast  oategory.  If  the  quantity  was  less  than 
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fig.  3.3.  500  mb  composite  height  differences  in  meters  of  Xaafc  minus 

alow  oy clones  in  the  Atlantic.  The  large  dot  marks 
the  position  of  the  cyclone.  (From  Xu  and  Gray,  1982). 


Fig.  3.4.  500  mb  composite  height  differences  in  meters  of  fast  minus 

slow  cyclones  in  the  Pacific.  The  large  dot  marks  the 
position  of  the  cyclone.  (From  Xu  and  Gray.  1982). 


zero,  the  long  wave  pattern  was  assumed  to  be  favorable  for  slow  motion 
so  the  cyclone  was  placed  in  the  slow  category.  It  was  hoped  that  this 
screening  would  be  able  to  omit  the  slow  moving  cyclones  (less  than  2.5 
m/s)  from  one  set  of  equations  and  the  fast  moving  cyclones  (greater 
than  7.5  m/s)  from  the  other  set  of  equations. 

In  a  conceptual  sense,  as  shown  in  Figs.  3.5  through  3.7,  different 
500  mb  long  wave  patterns  are  associated  with  fast  and  slow  motion  and 
it  was  felt  that  this  procedure  would  result  in  a  crude  classification 
according  to  synoptic  types  in  a  manner  similar  to  the  scheme  proposed 
by  Tse  (1966)  for  typhoon  forecasting  in  the  Pacific.  If  this  theory 
holds,  then  certain  predictors  would  be  selected  in  different  orders  or 
given  different  weights  for  the  fast  and  slow  cases. 

3.2  24-Hour  Time  Steps 

For  extended  forecast  periods,  a  likely  weakness  of  a 
stratification  based  on  the  position  of  the  cyclone  relative  to  the 
ridge  is  that  a  tropical  cylone's  position  relative  to  the  ridge  may 
change  over  a  span  of  three  days.  Thus,  although  the  initial  position 
of  the  cyclone  may  place  it,  for  example,  south  of  the  ridge,  it  may  be 
located  north  of  the  ridge  48  hours  later  and  the  environmental 
Influences  affecting  its  movement  over  the  next  24  hours  will  certainly 
differ  from  those  affecting  its  movement  had  it  remained  south  of  the 
ridge  during  the  entire  72-hour  period. 

In  order  to  eliminate  the  problem  of  a  cyclone  initially  in  one 
category  moving  into  another  category  during  the  oourae  of  a  72-hour 
motion,  and  also  to  attempt  to  better  describe  the  environment  in  which 
the  cyclone  was  embedded  as  it  changed  over  the  period  of  72  hours,  the 
72-hour  period  was  segmented  into  three  24-hour  periods.  The  segmenting 
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Fig.  3.5.  a)  Favorable  height  pattern  for  fast  motion  north  of  the 
subtropical  ridge,  b)  Typical  height  pattern  for  slow 
motion  north  of  the  subtropical  ridge.  (From  Xu  and 
Gray,  1982.) 
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Fig.  3.6.  a)  Typical  height  pattern  favorable  for  fast  notion  on  the 

subtropical  ridge,  b)  Typical  height  pattern  for  slow  motion 
on  the  subtropical  ridge.  (From  Xu  and  Gray,  1982.) 
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a)  Favorable  height  pattern  for  fast  motion  south  of  the 
subtropical  ridge,  b)  Typical  height  pattern  for  slow  motion 
south  of  the  subtropical  ridge.  (From  Xu  and  Gray,  1982.) 
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of  72-hour  notion  into  3  discrete  24-hourly  time  steps  departs  from  the 
traditional  statistical  approach.  There  is  an  important  distinction 
between  CSU84  and  the  other  statistical  models  in  forming  the 
developmental  set.  The  other  models  require  72-hour  tracks  for  each 
sample  case  whereas  each  sample  case  for  CSU84  is  composed  of  24-hour 
tracks.  Thus,  it  is  possible  to  get  many  more  cases  to  develop  the 
regression  equations  for  the  CSU84  model.  For  the  Atlantic,  there  are 
977  cases  for  the  south-of-ridge  category,  1034  cases  for  the  north-of- 
ridge  category  and  371  cases  for  the  on- th e-ridge  category.  For  the  NW 
Pacific,  there  are  4074  cases  south  of  the  ridge,  799  cases  north  of  the 
ridge  and  973  cases  on  the  ridge.  With  the  discretized  steps,  CSU84 
calculates  new  past-motion  persistence  predictors  and  updates  synoptic 
grid-point  data  relative  to  the  new  cyclone  position  at  the  start  of 
each  24-hour  forecast.  The  direction  of  motion  over  the  previous  24 
hours  or  a  subjective  analysis  of  the  new  position  of  the  cyclone 
relative  to  the  subtropical  ridge  is  used  to  determine  changeover  from 
one  stratification  class  to  another.  With  the  objective  method  of 
changing  categories  for  cyclones  in  the  south-of-ridge  category,  if  the 
motion  over  the  previous  24-hour  period  indicates  a  heading  between  330° 
and  30°,  that  segment  is  placed  in  the  on-the-ridge  category  for  the 
formulation  of  the  next  24-hour  motion.  Similarly,  for  cyclones  in  the 
south-of-ridge  or  on-the-ridge  categories,  if  the  motion  over  the 
previous  24-hour  period  is  between  31°  and  120°,  that  segment  is  placed 
in  the  north-of-ridge  category  for  the  formulation  of  the  next  24-hour 
motion.  The  perfect-prog  approach  is  taken  to  provide  the  new  synoptic 
information.  In  an  operational  environment,  the  24-  and  48-hour  500  mb 
height  and  surface  pressure  prognostic  fields  would  be  used  to  provide 
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the  information  to  formulate  the  48-hour  and  72-hour  forecast  motions, 
respectively.  In  other  words,  the  24-hour  prognostic  fields  are  applied 
to  compute  only  the  24  to  48  hour  motion  and  the  48-hour  prognostic 
fields  are  applied  to  compute  only  the  48  to  72  hour  motion.  This  is  a 
significant  departure  from  the  traditional  method  of  introducing 
prognostic  data.  Figure  3.8  is  a  flowchart  depicting  the  3-step 
procedure  for  arriving  at  the  24-,  48-  and  72-hour  forecast  positions. 

3.3  Application  of  Prognostic  Data 

The  method  of  applying  the  prognostic  data  in  CSU84  is  unique. 

Other  statistical  models  that  employ  prognostic  data  Invariably  use 
several  forecast  fields  (24-,  36-,  and  48-hour)  for  every  forecast 
period,  l.e.,  0  to  12,  0  to  24,  0  to  36,  0  to  48,  0  to  72.  The  weakness 
of  such  an  approach  is.  as  pointed  out  by  Neumann  (1975),  the  over¬ 
weighting  of  the  predictors  from  the  prognostic  data  fields  at  the 
expense  of  observed  data  fields.  This  is  understandable  when  the 
perfect-prog  approach  is  used  to  introduce  the  prognostic  data  since  all 
the  prog  data  would  be  "perfect”  —  the  data  come  from  verifying 
analysis  fields.  Thus,  it  is  of  no  surprise  that  these  models  degrade 
tremendously  when  actual  forecast  fields,  which  are  not  as  accurate  as 
the  perfect-prog  development  data,  are  used  to  introduce  the  synoptic 
predictors  in  operational  use.  This  point  will  be  amplified  in  Chapter  5. 

The  unique  method  of  applying  prognostic  data  developed  in  CSU84 
minimizes  this  degradation  in  two  ways.  First,  by  using  just  the 
current  and  24-hour  old  analyses  to  oompute  the  first  24  hours  of 
motion,  one  has  a  "solid"  24-hour  forecast  position  on  which  to  build 
a  48-hour  forecast  based  on  prognostic  data  to  arrive  at  a  72-hour 
forecast  position.  Second,  the  48-hour  forecast  position  is  computed 


Flow  Chart  To  Compute  The  24-HR,  48-HR.  And  72-HR  Forecast  Positions 


Fig.  3.8.  3-step  flow  chart  depicting  the  procedure  for  computing  the  24-,  48-  and  72-hour  forecast 
positions. 
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based  on  a  24-hour  motion  using  the  24-hour  forecast  fields  and  the 
initial  analysis  fields  to  provide  the  tendencies  of  heights  and 
pressures.  Hence,  just  one  forecast  field  is  required  to  compute  the 
second  24-hour  motion.  Thus,  it  is  only  the  third  and  final  24-hour 
motion  that  depends  on  two  forecast  fields  (24-hour  and  48-hour).  It 
would  appear  logical  to  shorten  the  forecast  Interval  that  a  scheme 
relies  on  as  much  as  possible  so  that  if  the  prog  field  is  in  error,  the 
error  in  the  ensuing  forecast  position  will  be  kept  to  a  minimum. 
Further,  although  one  or  both  of  the  forecast  fields  may  be  in  error, 
the  tendencies  may  still  be  fairly  accurate.  Thus,  by  discretizing  the 
72-hour  cyclone  motion  into  three  shorter  24-hour  segments,  the  reliance 
of  the  forecast  scheme  on  "perfect”  synoptic  data  to  predict  a  24-hour 
motion  should  not  be  as  unrealistic  as  the  reliance  it  would  have  if  it 
was  used  to  predict  a  72-hour  motion.  If  this  notion  holds,  then  the 
degradation  in  skill  one  can  expect  when  using  actual  prognostic  data 
should  not  be  as  drastic  with  this  scheme. 

This  CSU84  method  certainly  appears  to  be  superior  to  ones  that 
attempt  to  compute  a  24-hour  forecast  position  based  on  analyses  as  well 
as  prognostic  data  and  then  to  recompute  48  and  72  hours  of  motion 
always  from  the  initial  position  using  three  forecast  fields  (24-,  36-, 
and  48-hour),  all  of  which  oould  contain  significant  errors. 

A  schematic  diagram  of  the  three  steps  needed  for  a  forecast 
projection  of  72  hours  is  presented  in  Fig.  3.9.  Note  that  what  is 
meant  by  past  motion  and  current  and  24-hour  old  synoptic  data  changes 
as  each  24-hour  displacement  is  stepped  forward.  As  a  result,  the 
actual  values  of  the  selected  predictors  will  change  as  the  forecast 
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positions  are  projected  forward  in  time.  This  is  because  all  of  the 
predictors  are  relative  to  the  position  of  the  cyclone  at  the  time  they 
are  applied.  Hence,  prognostic  data  are  used  as  if  they  were  analysis 
data  at  their  respective  valid  times. 

The  cyclone  symbols  with  the  darkened  circle  are  the  positions  at 
which  the  three  24-hour  displacements  are  computed.  The  two  cyclone 
symbols  with  the  open  circle  to  the  left  of  these  positions  represent 
the  12-  and  24-hour  old  positions  relative  to  each  of  the  three 
positions  at  which  a  24-hour  forecast  is  made.  The  cyclone  symbol  with 
the  open  circle  to  the  right  of  the  darkened  circles  is  the  24-hour 
displacement.  As  the  72-hour  forecast  process  is  stepped  forward  in 
24-hour  Increments,  each  24-hour  displacement  position  becomes  the  new 
position  from  which  to  oompute  the  next  24-hour  displacement. 

The  "H"  and  "P"  above  the  cyclone  symbols  indicate  which  500  mb 
height  and  surface  pressure  fields  provide  the  required  synoptic 
predictors  to  compute  the  displacements.  With  Tq  the  time  at  the 
initial  position.  H(o)  and  P(o)  corresponds  to  the  current  analyses, 
H(-24)  and  P<-24)  correspond  to  24-hour-old  analyses,  H(+24)  and  P(+24) 
correspond  to  the  24-hour  prognoses  and  H<+48)  and  P(+48)  correspond  to 
the  48-hour  prognoses.  It  must  be  kept  in  mind  that  the  gridpolnt 
predictors  are  relative  to  the  position  of  the  cyclone  at  the  time  each 
forecast  is  made,  not  relative  to  the  initial  position  of  the  cyclone. 


4.  THE  REGRESSION  EQUATIONS 
4.1  The  Regression  Equation 

To  develop  the  prediction  equations,  a  stagewlse-stepwise  multiple 
linear  regression  procedure  proposed  by  Lund  (1971)  was  adopted.  Due  to 
the  large  number  of  both  sample  cases  and  potential  predictors,  the 
normal  stepwise  regression  approach  would  have  been  inappropriate  since 
just  the  two  covariance  arrays  containing  the  sums  of  squares  and 
cross-products  would  have  exceeded  the  memory  size  of  the  largest 
available  computers.  With  the  procedure  suggested  by  Lund,  there  is 
just  one  relative  coefficient  matrix  instead  of  two  covariance  matrices. 

The  regression  equations  were  obtained  using  various  computational 
subroutines  available  in  The  International  Mathematics  and  Statistics 
Libraries  (IMSL.  1982)  that  reside  on  the  NCAR  CRAY-1  computer.  The 
method  employed  was  as  follows: 

(1)  Means,  standard  deviations  and  correlation  coefficients  were 
computed  for  each  predictor 

(2)  IMSL  subroutine  BECOVM  computed  the  corrected  sum  of  squares 
and  cross-products 

(3)  IMSL  subroutine  RLMUL  performed  an  analysis  of  variance 
including  the  standard  deviation  of  the  residuals  and  the  regression 
coefficient  and  intercept  estimates. 

(4)  IMSL  subroutine  RLRES  performed  a  residual  analysis  on  the 
predictors  used  in  fitting  a  regression  model.  The  residual  from  this 


fitted  model  was  used  as  the  new  dependent  variable  (predictand)  and 
steps  (1)  -  (4)  repeated  20  times  to  select  the  20  best  predictors. 

(5)  IMSL  subroutine  RLSEP  performed  a  forward  stepwise  regression 
analysis  using  the  20  predictors  selected  from  the  stagewise  procedure 
to  compute  the  regression  coefficients.  The  process  that  is 
accomplished  by  the  first  four  steps  has  also  been  referred  to  as 
successive  elimination.  After  the  removal  of  each  selected  predictor,  a 
new  residual  is  calculated  with  the  remaining  predictors  and  that 
becomes  the  new  predictand  with  which  the  remaining  predictors  are 
correlated. 

4.2  Geographical  Inhomogeneity 

Northern  Hemisphere  cyclone  tracks  exhibit  significant  regional 
variations.  In  the  NW  Pacific,  the  near- equatorial  monsoon  trough  is 
responsible  for  easterly  or  erratic  motion  with  low  latitude  cyclones. 
When  cyclones  manage  to  break  through  the  subtropical  ridge,  they 
sometimes  travel  far  poleward  before  being  steered  eastward  by  the 
westerlies.  For  the  South  China  Sea  region,  complicated  and  conflicting 
wind  flow  patterns  between  the  southwest  monsoon  flow  and  the 
continental  northeaster lies  lead  to  a  strong  vertical  wind  shear  and 
extremely  erratic  motion  is  common.  In  the  North  Indian  Ocean,  weak 
steering  currents  lead  to  slow  motion.  In  the  Atlantic,  cyclones 
generally  track  in  a  west  to  northwestward  direction  recurving  at  higher 
latitudes.  The  ratio  of  very  fast  and  very  slow  moving  cyclones  to  the 
total  is  high,  however,  resulting  in  difficult  forecast  situations. 

Regression  equations  were  therefore  developed  separately  for  the 
four  geographical  regions  using  the  previous  set  of  potential  predictors 
but  with  slight  regional  changes.  Due  to  the  low  annual  frequency  of 


cyclones  in  the  North  Indian  Ocean,  it  was  necessary  to  use  500  mb 
computed  geostrophic  winds  in  place  of  200  mb  computed  geostrophic  winds 
since  200  mb  data  are  available  only  after  1962.  This  permitted  the  use 
of  the  entire  historical  best  tracks  going  back  to  1946.  In  addition, 
to  further  increase  the  sample  size,  no  speed  stratification  was 
attempted  in  the  North  Indian  Ocean. 

For  the  South  China  Sea,  since  the  vertical  shear  of  the  steering 
flow  appeared  to  be  of  overiding  importance  rather  than  the  location  of 
the  mid-level  subtropical  ridge  for  cyclone  motion  in  this  region,  these 
cyclones  were  not  stratified  relative  to  the  500  mb  subtropical  ridge. 
Storms  in  this  region  were  also  not  stratified  by  speed  both  because 
they  tend  to  move  slowly  and  also  in  order  to  obtain  as  large  a  sample 
size  as  possible.  Instead,  additional  steering  levels  -  700  mb  and  850 
mb  -  were  incorporated  in  an  attempt  to  better  depict  the  steering  flow 
at  different  vertical  levels. 

For  the  NW  Pacific,  two  observed  phenomena  presented  special 
difficulties  in  stratification.  One  was  the  erratic  motion  of  cyclones 
embedded  in  the  near-equatorial  monsoon  trough.  The  other  was  the  long 
northward  track  of  cyclones  breaking  through  the  subtropical  ridge  prior 
to  being  steered  eastward  by  the  upper  level  westerlies. 

In  the  former  case  of  anomalous  motion,  it  is  obvious  that  the 
direction  of  motion  is  a  poor  indicator  of  future  motion.  Instead  of 
deleting  these  monsoon-trough  cases,  these  cyclones  were  left  intact  in 
the  developmental  data  set  with  the  understanding  that  correlation  of 
motion  south  of  the  ridge  with  synoptic  predictors  would  not  be  as  high 
as  those  for  the  Atlantic  cyclones.  To  handle  the  long  northward  trek 
of  cyclones  in  the  Pacific,  the  criteria  for  the  on-the-ridge  category 
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was  broadened  to  include  all  segments  that  had  a  direction  of  motion 
between  330°  and  30°  rather  than  restricting  this  category  to  24  hours 
before  and  after  the  recurvature  point. 

4.3  The  Potential  Predictors 

Equations  for  the  east-west  and  north-south  components  of  cyclone 
motion  were  developed  by  screening  over  700  potential  predictors.  The 
potential  predictors  included  the  past  12-hour  and  past  24-hour  east- 
west  and  north-south  motions  and  grid  point  synoptic  data  of  "current" 
and  the  past  24-hour  changes  of  S00  mb  heights  and  sea  level  pressures. 
Recall  that  all  of  the  predictors  are  applied  relative  to  a  new  cyclone 
position  each  time  a  24-hour  motion  is  computed.  Thus,  the  synoptic 
grid  and  the  12-  and  24-hour  old  cyclone  positions  needed  to  calculate 
the  past  motions  move  with  the  cyclone. 

As  a  result,  the  analysis  field  valid  at  the  initial  time  Tq  is  the 
"current"  synoptic  field  and  that  field,  along  with  the  analysis  field 
valid  at  Tq  minus  24  hours  (T_24> .  are  used  to  calculate  the  24-hour 
changes  of  S00  mb  heights  and  sea  level  pressures.  These  two  fields  are 
also  used  to  compute  the  "current"  and  24-hour  changes  of  the  steering 
predictors.  Once  the  first  24-hour  displacement  has  been  stepped 
forward,  the  24-hour  forecast  field  valid  at  Tq  plus  24  hours  (T+24) 
becomes  the  "current"  synoptic  field  and  the  analysis  field  at  Tq 
becomes  the  24-hour-old  field  and  all  of  the  grid-point  predictors  as 
well  as  the  updated  persistence  predictors  are  computed  relative  to  the 
new  cyclone  position  at  T+2^.  When  the  second  24-hour  displacement  is 
stepped  forward  from  the  T+2^  cyclone  position,  the  48-hour  forecast 
field  valid  at  Tq  plus  48  hours  (T+4g)  becomes  the  "current"  synoptic 
field  and  the  forecast  field  valid  at  T+24  becomes  the  24-hour-old 
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field.  Then,  all  of  the  grid-point  predictors  as  well  as  the  updated 
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persistence  predictors  are  computed  relative  to  the  new  cyclone  position 
at  T+4g.  Using  these  updated  predictors,  the  final  24-hour  displacement 
is  stepped  forward  from  the  T+4g  cyclone  position  to  determine  the  72- 
hour  forecast  position.  All  references  to  steering  winds  apply  to 
geostrophic  winds  which  were  calculated  from  the  gradients  of 
geopotential  heights  using  the  following  formula: 

For  the  wind  at  the  cyclone  center,  the  U  and  V  components  are 
computed  as: 


U  “  (<H-10,~10  +  H-10.-5  +  H-10,0  +  H-10, 5  +  H-10.10) 

_(H10,-10  +  H10,-5  +  H10,0  +  H10, 5  +  *10.10* *  / 

(5  x  sin  p) 


V  =  <<H-10,10  +  H-5,10  +  H0.10  +  H5,10  +  H10.10) 
"(H-10.-10  +  H-5,-10  +  H0,-10  +  H5,-10  +  H10,-10))  1 


(5  x  sin  (p)  x  cos  (p)) 

For  the  wind  10°  northwest  of  the  cyclone  center,  the  U  and  V 
components  are  computed  as: 


U  ■  <H5.-10  -  “l5.-10>  '  •*"  <*> 

v  ■  "  H10,-15*  ^  ^Sin  W  x  003  ^ 

For  the  wind  5°  northwest  of  the  cyclone  center,  the  U  and  V 

components  are  computed  as: 

U  ”  <H0,-5  -  H10,-5>  '  sin  <*> 

V  ■  (H5#o  "  H5,-10)  ^  <sin  ^  x  008  W* 

For  the  wind  10°  north  of  the  cyclone  center,  the  U  and  V 

components  are  computed  as: 

0  “  <H5,0  '  H15,0)  1  8in  * 

V  “  <H10  S  ~  H10  -«>  t  (sin  *  x  008  & 


For  the  wind  5  north  of  the  cyclone  center,  the  U  and  V  components 


are  computed  as: 

U  =  Q  -  H10  q)  /  sin  ? 

V  =  (H5  s  -  h5  _5)  /  (sin  0  x  cos  £) 

For  the  wind  15°  northwest  of  the  cyclone  center,  the  U  and  V 
components  are  computed  as: 

U  =  (H10,-15  "  H20,-15)  /  sin  * 

V  =  (H15.-10  "  H15,-20)  1  (sin  *  x  003  9) 

For  the  wind  20°  northwest  of  the  cyclone  center,  the  U  and  V 

components  are  computed  as: 

U  =  (H15,-20  -  H25,-20)  '  3in  * 

V  “  (H20,-15  "  H20,-25)  1  (3in  *  x  003  0> 

For  the  wind  15°  north  of  the  cyclone  center,  the  U  and  V 

components  are  computed  as: 

U  =  (H10.0  "  H20,0}  /  3in  9 

V  =  ^H15,5  "  H15  -5*  ^  <sin  9  x  003  9^ 

For  the  wind  20°  north  of  the  cyclone  center,  the  U  and  V 

components  are  computed  as: 

U  =  (H15.0  "  H25.0)  /  3in  9 

V  =  (H20.5  "  H20,-5>  '  <3in  9  x  003  9) 

j  is  the  S00  mb  or  200  mb  height  at  grid  point  i,j  where  i 
represents  latitude  and  J  represents  longitude  from  the  cyclone  center. 
Latitudes  are  positive  to  the  north  and  longitudes  are  positive  to  the 
east  of  the  cyclone  center,  p  is  the  latitude  of  the  grid  point  at 
which  the  steering  wind  is  being  computed. 


Table  4.1  lists  the  732  potential  predictors  for  the  Atlantic  and 
NW  Pacific.  Table  4.2  lists  the  728  potential  predictors  for  the  North 
Indian  Ocean  and  Table  4.3  lists  the  1364  potential  preditors  for  the 
South  China  Sea.  The  list  of  728  potential  predictors  for  the  Atlantic 
and  NW  Pacific  using  500  mb  steering  instead  of  200  mb  steering  is 
identical  to  those  listed  for  the  North  Indian  Ocean.  The  500  mb 
steering  wind  over  the  cyclone  center  is  computed  for  all  versions  of 
the  model.  This  is  the  reason  for  the  difference  in  the  number  of 
predictors  between  the  two  versions  for  the  NW  Pacific  and  the  Atlantic. 

Variations  of  the  model  using  surface  pressure  gradients  as  well  as 
thickness  of  the  surface  to  500  mb  surfaces  were  tested.  However,  the 
screening  program  failed  to  select  either  category  of  predictors. 

The  three  categories  of  predictors  used  in  the  CSU84  prediction 
model  were  shown  in  Chapter  3.  As  explained  earlier,  the  "current" 
and  the  24-hour-old  synoptic  data  are  not  always  the  analysis  and  the 
24-hour-old  analysis  at  the  initial  time  for  the  24-  to  48-  and  48-  to 
72-hour  predictions.  These  fields  change  as  the  three  24-hour 
displacements  are  stepped  forward.  In  development,  a  perfect-prog 
approach  was  used  to  introduce  analysis  data  valid  at  the  verifying 
times.  In  actual  practice,  24-  and  48-hour  forecast  fields  would  be 
used  at  their  verifying  times. 

Importance  of  Different  Predictors.  The  reduction  of  variance  due 
to  each  of  the  three  categories  of  predictors  was  examined  for  Atlantic 
and  Paific.  A  few  surprising  observations  are  evident  from  Table  4.4. 
One  was  the  significant  difference  in  the  variance  reducing  potential  of 
the  steering  component  for  the  north-of-ridge  and  on-the-ridge  between 


the  Atlantic  and  NW  Pacific 


TABLE  4.1 


The  732  potential  forecast  predictors  for  the  Atlantic  and  NW  Pacific 
1-170  24-hr  500  mb  height  changes  at  the  17  X  10  grid  points  [AH] 

171-340  500  mb  heights  at  the  17  X  10  grid  points  [H] 

341-510  24-hr  sea  level  pressure  changes  at  the  17  X  10  grid  points 

(mb)  [AP] 

511-680  Sea  level  pressures  at  the  17  X  10  grid  points  (mb)  [P] 

681-682  Past  24-hr  0  and  V  cyclone  motion  (km/hr)  [U^,  V_24] 

683-684  Past  12-hr  0  and  V  cyclone  motion  (km/hr)  [U_12,  V_12l 

685-686  500  mb  0  and  V  wind  components  at  cyclone  center  If500u*  f500y] 

687-688  24-hr  500  mb  U  and  V  wind  change  at  cyclone  center  [Af500u. 

^OOv1 

689-690  200  mb  0  and  T  wind  components  at  cyclone  center  t^oou'  $200v^ 

691-692  24-hr  200  mb  U  and  T  wind  change  at  cyclone  center  [Af200u* 

A*200vI 

693-694  200  mb  0  and  V  wind  components  10°  northwest  of  cyclone 

center  If2ooNW10u’  ^OONWlOv1 

695-696  24-hr  200  mb  0  and  V  wind  change  10°  northwest  of  cyclone 

center  [At2ooNW10u*  A*200NW10v] 

697-698  200  mb  U  and  V  wind  components  5°  northwest  of  cyclone 

center  lf2ooNW5u'  *200NW5v] 

699-700  24-hr  200  mb  0  and  V  wind  change  5°  northwest  of  cyclone 

center  lAf200NW5u»  A*200HW5v1 

701-702  200  mb  0  and  V  wind  difference  between  10°  northwest  and 

5°  northwest  of  cyclone  center  [»200Nwio-5u'  *200NV10-5v1 

703-704  24-hr  200  mb  U  and  V  wind  change  difference  between  10° 

northwest  and  5°  northwest  of  cyclone  center  lAf200NW10-5u' 

A*200NW10-5v1 

705-706  200  mb  U  and  V  wind  components  10°  north  of  cyclone 

center  [*2ooNlOu*  *200N10v1 

707-708  24-hr  200  mb  U  and  V  wind  change  10°  north  of  oydone 

center  IAf2ooNlOu'  A*200N10vJ 
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709-710 

711-712 

713-714 

715-716 

717-718 

719-720 

721-722 

723-724 

725-726 

727-728 

729-730 

731-732 


TABLE  4.1  (cont'd) 

200  mb  U  and  V  wind  components  5°  north  of  cyclone  center 
I#200N5u'  1200N5V1 

24-hr  200  mb  U  and  V  wind  change  5°  north  of  cyclone  center 
IA*200N5u'  A*200N5v] 

200  mb  U  and  V  wind  difference  between  10°  north  and  5° 
north  of  cyclone  center  [f200N10_5u.  ^ooNlO-Sv1 

24-hr  200  mb  U  and  V  wind  change  difference  between  10°  north 
and  5°  north  of  cyclone  center  tA?200N10_5u.  A} 200N10_5v 1 

200  mb  U  and  V  wind  components  15°  northwest  of  cyclone 
center  tf200NW15u.  *200NW15v] 

24-hr  200  mb  U  and  V  wind  change  15°  northwest  of  cyclone 
centers  I4*20omil5u.  MjOOKUlSv1 

200  mb  U  and  V  wind  components  20°  northwest  of  cyclone 
centers  If200NW20u*  *200NW20v] 

24-hr  200  mb  0  and  V  wind  change  20°  northwest  of  cyclone 
center  IA}200NW20u,  A#200NW20vl 

200  mb  U  and  V  wind  components  15°  north  of  cyclone  center 
I*200N15u'  1200N15V1 

24-hr  200  mb  U  and  V  wind  change  15°  north  of  cyclone  center 
tA*200N15u*  A*200N15vI 

200  mb  U  and  V  wind  components  20°  north  of  cyclone  center 
I#200N20u'  *200N20v3 

24-hr  200  mb  U  and  V  wind  change  20°  north  of  cyclone  center 
IA*200N20u'  A*200N20vJ 
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TABLE  4.2 

The  728  Potential  Forecast  Predictors  For  The  North  Indian  Ocean 
1-170  24-hr  500  mb  height  change  at  the  (17  X  10)  grid  points  [AH] 

171-340  500  mb  heights  at  the  17  X  10  grid  points  [H] 

341-510  24-hr  sea  level  pressure  change  at  the  17  X  10  grid  points  (mb) 

[AP] 

511-680  Sea  level  pressures  at  the  17  X  10  grid  points  (mb)  [P] 

681-682  Past  24-hr  0  and  V  cyclone  motion  (km/hr)  [U_2<>  V_24l 

683-684  Past  12-hr  0  and  V  cyclone  motion  (km/hr)  [U_12.  V_12l 

685-686  500  mb  U  and  V  wind  components  at  cyclone  center  lt5Q0u' 

•500V1 

687-688  24-hr  500  mb  U  and  V  wind  change  at  cyclone  center  U*500u' 

A*500vJ 

689-690  500  mb  0  and  T  wind  components  10°  northwest  of  cyclone 

center  If500NW10u.  fsoONWIOv1 

691-692  24-hr  500  mb  U  and  V  wind  change  10°  northwest  of  cyclone 

center  Uf5001wl0u.  At5ooNW10v] 

693-694  500  mb  0  and  V  wind  components  5°  northwest  of  cyclone 

center  [fS00NW5u.  fsoONWSv1 

695-696  24-hr  500  mb  U  and  V  wind  change  5°  northwest  of  cyclone 

center  IA1500NW5u.  ^fsoONWSv1 

697-698  500  mb  U  and  V  wind  difference  between  10°  northwest  and 

5°  northwest  of  cyclone  center  If500Nwl0_5u.  fsooNWlfr-Sv1 

699-700  24-hr  500  mb  U  and  V  wind  change  diffeence  between  10° 

northwest  and  5°  northwest  of  cyclone  center  IAl500NV10-5u’ 

^OONWIO-Sv1 

701-702  500  mb  U  and  V  wind  components  10°  north  of  cyclone  center 

^SOONlOu*  •500N10V1 

703-704  24-hr  500  mb  (J  and  V  wind  change  10°  north  of  cyclone  center 

lA*500N10u'  A*500N10v1 

705-706  500  mb  U  and  V  wind  components  5°  north  of  cyclone  center 

^•500N5u’  1500N5V1 


TABLE  4.2  (cont'd) 


707-708 

709-710 

711-712 

713-714 

715-716 

717-718 

719-720 

721-722 

723-724 

725-726 

727-728 


24-hr  500  mb  U  and  V  wind  change  5°  north  of  cyclone  center 
[A*500N5u'  A*500N5v] 

500  mb  U  and  V  wind  difference  between  10°  north  and  5° 
north  of  cyclone  center  lf500N10_5u.  fsooNio-Sv1 


24-hr  500  mb  U  and  V  wind  change  difference  between  10 

and  5°  north  of  cyclone  center  [A® ,Amj1 

*500N10-5u,  Af50()N10_5v 


north 

] 


500  mb  U  and  V  wind  components  15°  northwest  of  cyclone 
center  [f500NW15u-  f500NW15v] 


24-hr  500  mb  U  and  V  wind  change  15°  northwest  of  cyclone 
center  lAf500m5u.  Af500NW15v] 


500  mb  U  and  V  wind  components  20°  northwest  of  cyclone 
center  [f500NW20u'  *500NW20v1 


24-hr  500  mb  U  and  V  wind  change  20°  northwest  of  cyclone 
center  [Af5()0NW20u,  Af500NW20v] 

500  mb  U  and  V  wind  components  15°  north  of  cyclone  center 
If500N15u*  ^500N15v^ 

24-hr  500  mb  U  and  V  wind  change  15°  north  of  cyclone  center 
lAf500N15u*  A*500N15v] 

500  mb  U  and  V  wind  components  20°  north  of  cyclone  center 
[*500N20u’  *500N20v1 

24-hr  500  mb  0  and  V  wind  change  20°  north  of  cyclone 
center  [Af500N2()u,  Af5Q0N20v] 


TABLE  4.3 


The  Potential  Forecast  Predictors  For  The  South  China  Sea 

1-170  24-hr  500  mb  height  changes  at  the  17  X  10  grid  points  [AH] 

171-340  500  mb  heights  at  the  17  X  10  grid  points  [H] 

341-510  24-hr  sea  level  pressure  changes  at  the  17  X  10  grid  points 

(mb)  [AP] 


511-680 

681-850 

851-1020 

1021-1190 

1191-1360 

1361-1362 

1363-1364 


Sea  level  pressures  at  the  17  X  10  grid  points  (mb)  [P] 

24-hr  700  mb  height  changes  at  the  17  X  10  grid  points  [AH7] 

700  mb  heights  at  the  17  X  10  grid  points  [H7] 

24-hr  850  mb  height  changes  at  the  17  X  10  grid  points  [AH8] 

850  mb  heights  at  the  17  X  10  grid  points  [H8] 

Past  24-hr  0  and  V  cyclone  motion  (km/hr)  [0  ...  V 

“24  —24 

Past  12-hr  U  and  V  cyclone  notion  (km/hr)  [U_12# 
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TABLE  4.4 


Percent  of  variance  reduction  due  to  the  three  categories  of  predictors 
(see  Fig.  3.7)  for  the  Atlantic  and  the  NW  Pacific.  Version  CSU84-Q0  o: 
the  model  uses  500  mb  steering.  Version  CSU842Q0  uses  200  mb  steering. 
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For  the  Atlantic,  the  variance  reduction  due  to  steering  was 
overwhelming  for  the  meridional  motion  (V)  nortn  of  the  ridge  and  for 
both  components  of  motion  on  the  ridge.  For  the  NW  Pacific,  however, 
steering  contributed  little  to  the  reduction  of  variance,  with  one 
exception.  Nearly  50%  of  the  variance  of  zonal  motion  was  explained  by 
200  mb  steering  predictors  north  of  the  ridge. 

A  second  surprise  is  the  amount  of  variance  reduction  due  to 
synoptic  predictors  as  compared  to  persistence.  With  the  exception  of 
the  south-of-ridge  category,  the  synoptic  predictors  are  responsible  for 
as  much,  if  not  more,  of  the  variance  reduction  as  is  persistence  in 
both  ocean  basins.  This  certainly  indicates  the  model's  potential  for 
predicting  anomalous  cyclone  motion  if  an  anomaly  is  present  in  the 
analysis  or  forecast  fields. 

4.4  Correlation  Coefficient  Fields 

Figures  4.1  through  4.4  present  the  correlation  coefficient  fields 
between  the  500  mb  geopotential  heights/ sea  level  pressures  and  tropical 
cyclone  motion  in  the  Atlantic  north  of  the  ridge  prior  to  the  selection 
and  removal  of  any  predictors.  Northward  and  eastward  cyclone  motions 
are  positive.  The  correlation  coefficient  fields  for  the  other  two 
categories  relative  to  the  ridge  are  presented  in  Appendix  B. 
Correlations  between  cyclone  motion  and  the  synoptic  fields  are  stronger 
in  the  Atlantic  basin  than  in  any  of  the  other  basins. 

For  the  south-of-ridge  category  the  difference  in  the  degree  of 
correlation  between  the  synoptic  fields  and  cyclone  motion  for  the 
Atlantic  as  oompared  to  the  Pacific  is  particularly  striking.  Again, 
the  erratic  motion  of  Pacific  cyclones  embedded  in  the  monsoonal  trough 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  north-or-ridge  fast  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (o)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressur 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( f  )  and  the  stars  Indicate  the  gridpolnts  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  north-of-ridge  fast  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(  # )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  north-of-ridge  slow  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(#  )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  north-of -ridge  slow  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( S  )  and  the  stars  indicate  the  grldpolnts  at  which  predictors  were  selected. 
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appears  to  account  for  the  poor  correlation  of  cyclone  motion  with  the 
synoptic  fields  in  the  Pacific,  especially  south  of  the  ridge. 

In  general,  for  the  south-of-ridge  category,  above  normal  heights 
and  pressures  to  the  north  are  associated  with  westward  motion  (Figs. 

B.5  and  B.7,  panels  a  and  c),  and  above  normal  heights  and  pressures  to 
the  east  are  associated  with  northward  motion  (Figs.  B.6  and  B.8,  panels 
a  and  c).  The  24-hour  change  or  tendency  fields  (panels  b  and  d), 
although  weaker,  imply  the  same  correlation. 

For  the  north-of-ridge  category,  there  is  again  nothing  unexpected 
in  the  figures.  Eastward  motion  is  associated  with  falling  heights  and 
pressures  to  the  north  (Figs.  4.1  and  4.3)  and  northward  motion  is 
associated  with  falling  heights  and  pressures  to  the  north  and  west 
(Figs.  4.2  and  4.4).  The  correlations  are  more  concentrated  with 
stronger  gradients  of  correlation  than  for  the  south-of-ridge 
categories. 

The  ability  of  the  screening  process  to  pick  out  the  subtle 
differences  in  the  longitudinal  position  of  the  trough  between  the  fast 
and  slow  cases  is  encouraging.  Note  that  the  centers  of  maximum 
correlation  indicate  that  for  the  slow  cases,  the  trough  has  already 
passed  the  cyclone  and  is  located  east  of  the  cyclone.  Compare  the 
centers  of  maximum  correlation  relative  to  the  cyclone  for  the  four 
panels  between  Fig.  4.1  and  4.3. 

For  the  on-the-ridge,  the  correlation  fields  are  slightly  more 
complex.  In  general,  falling  heights  and  pressures  to  the  north  and 
rising  heights  and  pressures  to  the  south  are  associated  with  eastward 
motion  (Figs.  B.l  and  B.3);  falling  heights  and  pressures  to  the 


west  and  rising  heights  and  pressures  to  the  east  are  associated  with 
northward  motion  (Figs.  B.2  and  B.4). 

Figures  4.5  through  4.8  present  the  correlation  coefficient  fields 
between  the  500  mb  geopotential  heights/sea  level  pressures  and  tropical 
cyclone  motion  in  the  NW  Pacific  north  of  the  ridge  prior  to  the 
selection  and  removal  of  any  predictors.  The  correlation  coefficient 
fields  for  the  other  two  categories  relative  to  the  ridge  are  presented 
in  Appendix  C.  The  relatively  weak  correlation  between  the  synoptic 
fields  and  cyclone  motion  south  of  the  ridge  (Figs.  C.5  through  C.8)  is 
to  be  noted.  Most  of  the  variance  reducing  potential  comes  from  the 
12-hour  past  motion.  Synoptic  information  reduces  the  variance  by  only 
an  additional  10%.  This  explains  why  statistical  techniques  such  as 
HPAC  that  uses  only  climatology  and  persistence  perform  well  in  the  low 
latitudes.  The  surprisingly  weak  correlation  of  cyclone  motion  to 
synoptic  fields  in  this  category  for  the  NW  Pacific  can  be  explained  by 
the  presence  of  the  near-equatorial  monsoon  trough  in  this  part  of  the 
globe.  When  cyclones  develop  within  this  trough  region,  erratic  motion 
is  observed  as  the  cyclones  drift  aimlessly  due  to  the  lack  of  any 
definite  steering  flow  until  they  drift  poleward  and  are  picked  up  by 
the  westerlies. 

In  general,  correlations  are  similar  to  the  Atlantic  although  there 
are  differences.  One  of  the  most  noticeable  differences  is  the  sea 
level  pressure  fields  for  the  north-of-ridge  category.  Unlike  the 
Atlantic,  the  strong  correlation  between  low  pressures  to  the  north  and 
eastward  motion  is  missing.  Refer  to  panel  c  on  Figs.  4.5  and  4.7. 

Also  missing  is  the  strong  correlation  between  low  pressures  to  the 
northwest  and  northward  motion.  Refer  to  panel  c  on  Figs.  4.6  and  4.8. 
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Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  north-of-ridge  fast  0  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(f  )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Fig.  4.6.  Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  north-of-ridge  fast  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( #  )  and  the  stars  Indicate  the  gridpolnts  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  north-of-ridge  slow  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( $  )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  north-of-ridge  slow  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( #  )  and  the  stars  indicate  the  gridpolnts  at  which  predictors  were  selected. 
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The  strongest  correlation  occurs  between  the  500  mb  heights  and  fast 
zonal  motion  for  both  north-of-ridge  (Fig.  4.5  panel  a)  and  on-the-ridge 
categories  (Fig.  C.l  panel  a).  Note  that  although  the  center  of  the 
strongest  negative  correlation  10°  north -northwest  of  the  cyclone  center 
was  selected  as  one  of  the  predictors  for  fast  zonal  motion  for  the  on- 
the-ridge  category  (Fig.  C.l  panel  a),  the  center  of  the  strongest 
negative  correlation  10°  north  of  the  cyclone  center  for  the  north-of- 
ridge  category  (Fig.  4.5  panel  a)  was  not  selected.  The  difference  in 
the  selection/non-seleotion  is  attributable  to  the  order  in  which  the 
predictors  are  selected.  Recall  that  the  correlation  fields  shown  are 
the  correlation  prior  to  the  selection  and  removal  of  any  predictors. 

For  the  on-the-ridge  category,  the  500  mb  height  at  the  center  of  the 
strongest  correlation  was  selected  first  while  for  the  north-of-ridge 
category,  other  predictors  such  as  steering  and  24-hour  height  change 
were  selected  first.  When  this  occurred,  the  residual  that  remained  did 
not  correlate  well  with  the  original  height  correlation  field. 

North  .Jn.d,i»n  pgean 

Figures  4.9  and  4.10  present  the  correlation  coefficient  fields 
between  the  500  mb  geopotential  heights/sea  level  pressures  and  tropical 
cyclone  motion  in  the  North  Indian  Ocean  north  of  the  ridge.  The 
correlation  fields  for  the  other  two  categories  relative  to  the  ridge 
are  presented  in  Appendix  D.  Correlations  are  not  nearly  as  strong  as 
in  the  other  ocean  basins,  particularly  south  of  the  subtropical  ridge 
as  depicted  by  Figs.  D.3  and  D.4.  Neither  the  500  mb  heights  nor  the 
sea  level  pressure  correlate  very  strongly  with  cyclone  motion  south  of 
the  ridge.  This  explains  the  rather  poor  reduction  in  variance  —  only 
55%  for  zonal  motion  and  47%  for  meridional  motion  —  for  this  category 
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Pig.  4.9.  Correlation  coefficient  fields  (x  .01)  for  the  North  Indian  Ocean  nor th-of- ridge  0  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(4)  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 


Correlation  coefficient  fields  (x  .01)  for  the  North  Indian  Ocean  north-of-ridge  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
( 4 )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 


of  cyclones.  Although  the  reduction  in  variance  is  relatively  low,  the 
forecast  errors  for  the  cyclones  south  of  the  ridge,  as  will  be  shown  in 
the  next  chapter,  are  surprisingly  small  for  both  the  dependent  and 
independent  set  due  to  the  small  standard  deviation.  This  is  a  result 
of  the  generally  slow  movement  of  these  cyclone  systems. 

Strongest  correlation  for  this  ocean  basin  occurs  between  the  500 
mb  heights  and  zoned  and  meridional  motion  north  of  the  ridge  (panel  a 
on  Figs.  4.9  and  4.10).  Again,  the  reason  for  the  selection  of  the 
grid-point  nearest  the  heights  20°  east-northeast  of  the  cyclone  center 
that  exhibited  the  strongest  positive  correlation  with  meridional  motion 
north  of  the  ridge  and  the  non-selection  of  the  grid-point  nearest  the 
heights  10°  north  of  the  cyclone  center  that  exhibited  strong  negative 
correlation  with  zonal  motion  is  due  to  the  order  of  predictor 
selection. 

In  general,  grid  points  near  the  centers  of  strongest  positive  or 
negative  correlation  were  selected  to  provide  height  and  pressure 
information  to  empirically  predict  24-hour  cyclone  motion.  Figure  4.10 
shows  a  good  example  of  the  correlation  field  between  24-hour  sea  level 
pressure  change  and  meridional  motion  north  of  the  ridge.  Refer  to 
panel  d  on  Fig.  4.10. 

Table  4.5  illustrates  the  reduction  of  variance  by  the  various 
predictors  for  the  North  Indian  Ocean.  Unlike  the  Atlantic  and  NW 
Pacific,  steering  does  not  contribute  to  cyclone  motion  except  for  a 
slight  contribution  to  the  zonal  motion  north  of  the  ridge.  Instead,  it 
is  the  synoptic  predictors,  primarily  the  surface  pressure  tendencies, 
that  play  a  significant  role.  With  the  exception  of  the  south-of-ridge 
category,  th«y  are  dominant.  Such  a  large  difference  from  the  NW 
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TABLE  4.5 

Percent  of  variance  reduction  by  the  various  predictors  for  the  North 
Indian  Ocean.  (X  -  persistence;  H  -  500  mb  height;  AH  -  24-hour  change 
of  500  mb  height;  P  -  surface  pressure;  AP  -  24-hour  change  of  surface 
pressure;  S  -  steering  at  500  mb). 

12  1 


X  -  31 
H  «  5 

North  AH  -  4  42 

(237  cases)  P  ■  6 

AP  -  14 
S  -  14 


X  -  18 
H  *  17 
AH  -  7  56 

P  -  4 
AP  «  28 
S  =  0 


X  -  29 
H  -  5 

On  AH  -  6  46 

(188  cases)  P  *  14 
AP  -  21 
S  *  5 


X  -  32 
H  =  0 

AH  *  9  27 

P  -  0 

AP  <*  18 
S  «  6 


X  -  41 

H  *=  7  | 

South  AH  -  2  15 

(188  cases)  P  *  1  ( 

AP  «  5  ' 

S  -  1 


X  «  33 
H  -  5 

AH  =  1  13 

P  =  3 
AP  «  4 
S  -  1 


Pacific  may  explain  why  an  attaspt  to  use  regression  equations  developed 
for  NW  Pacific  cyclones  failed  when  used  for  North  Indian  Ocean 
cyclones. 

The  screening  process  was  terminated  when  a  selected  predictor 
contributed  less  than  1  percent  to  the  reduction  in  the  variance  or  when 
a  maximum  number  of  20  predictors  had  been  selected.  As  will  be 
explained  in  Chapter  5.  speed  stratification  failed  to  enhance  the 
predictive  skill  of  the  model.  As  a  result,  the  regression  equations 
presented  in  this  chapter  apply  to  the  model  versions  with  no  speed 


stratification. 
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Tables  4.6  through  4.8  list  the  selected  predictors,  the  associated 
coefficients  and  the  reduction  of  variance  for  the  three  stratified 
categories  in  the  Atlantic.  Similar  information  for  the  NW  Pacific  is 
presented  in  Tables  4.9  through  4.11.  Tables  4.12  through  4.14  contain 
the  list  for  the  North  Indian  Ocean  and  Table  4.15  for  the  South  China 
Sea. 

Note  that  for  all  of  the  ocean  basins,  the  majority  of  the  selected 
synoptic  predictors  are  located  at  grid  points  15°  or  greater  from  the 
cyclone  center.  Many  of  them  are  located  on  the  edge  of  the  grid  40° 
from  the  cyclone.  It  was  gratifying  to  observe  that  the  large-scale 
circulation  features  were  being  incorporated  into  the  model  without 
explicitly  forcing  out  of  the  regression  any  grid  points  near  the 
cyclone. 

Further,  as  many,  if  not  more,  24-hour  past  changes  of  the  500  mb 
heights  and  sea  level  pressure  were  selected  compared  to  the  current 
heights  or  pressures.  Apparently,  by  not  using  prognostic  data  prior  to 
its  valid  time,  the  tendencies  of  heights  and  pressures  were  critical  in 
determining  how  the  current  fields  would  evolve  over  the  next  24  hours. 
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5.  VERIFICATION  AND  RESULTS 


5.1  Atlantic 

For  easy  reference.  Table  5.1  presents  a  complete  list  of  the 
various  model  runs  that  were  performed  with  Atlantic  cyclones.  The 
symbolic  name  for  each  run  was  selected  to  easily  denote  the  key 
elements  for  each  particular  model  run.  For  example,  the  first  run  is 
deciphered  as: 


500mb 
steering 


CSU84  ♦  / 

X 


dependent  set 

no  progs 


,speed  stratification 


C5  -  D  -  NP  (s)- 
46-,78/BTs. 


1946 
data 


Y  \ 


best  track 
initial  position 


For  each  run,  one  or  more  of  the  following  elements  were  changed: 

1)  500  mb  or  200  mb  steering  level  (5  or  2) 

2)  dependent  or  Independent  data  set  or  both  combined 

3)  perfect-prog  (PP)  (verifying  analyses),  actual  progs  (AP) 

or  no  progs  (NP) 

4)  speed  stratification  (s)  or  no  speed  stratification  (blank) 

5)  data  period  of  1946-78,  1962-81,  1968-81,  1973-79,  or  1979-81 

6)  best  track  (BT)  Initial  positions  or  perturbed  best  track  (PT) 

initial  positions. 

All  of  the  forecast  displacement  errors  referred  to  in  this  paper  are  in 
nautical  miles  (n  ml)  and  this  mean  vector  error  (great-circle  distance 
between  the  forecast  and  the  actual  location)  was  computed  using  the 
following  formula: 

E  -  cos'1  (  sin  XLAT  x  sin  FLAT  +  cos  XL AT  x  cos  FLAT 
x  cos  (XL ON  -  FLON)  ]  x  60 
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where  XLAT  and  XLON  are  the  best  track  latitude  and  longitude  and  PLAT 
and  FLON  are  the  forecast  latitude  and  longitude.  The  result  of  the 
verification  on  the  dependent  data  set  used  to  develop  the  regressions 
for  the  Atlantic  is  shown  in  Tables  5.2  and  5.3.  The  three  model  runs 
shown  in  Table  5.2  used  500  mb  steering  predictors  instead  of  200  mb 
steering  predictors  with  the  cases  categorized  with  respect  to  the 
subtropical  ridge  and  further  stratified  into  fast  and  slow  subsets 
based  on  the  500  mb  height  patterns.  C5-D-NP(s)/ (46-78/BT)  is  the 
traditional  or  classical  approach  of  applying  different  regression 
equations  for  each  of  the  three  cyclone  displacements  —  0  to  24.  0  to 
48.  and  0  to  72.  For  the  synoptic  predictors,  each  of  the  equations  use 
only  the  current  and  24-hour  old  analyses  with  no  prognostic  data.  C5- 
D-PP(s) /(46-78/BT)  used  the  same  data  set  as  C5-D-NP(s) / (46-78/BT)  and 
introduced  the  new  approach  of  segmenting  the  forecast  motion  into  three 
24-hour  steps.  It  applied  the  same  regression  equations  three  times  as 
the  forecast  positions  were  stepped  forward.  It  used  perfect  prog  data 
to  update  the  synoptic  predictors  relative  to  the  new  cyclone  position 
at  the  24-hour  and  48-hour  forecast  positions.  C5-D-PP(s)/(62-81/BT) 
applied  the  same  technique  as  C5-D-PP(s)/(46-78/BT)  to  a  smaller,  more 
recent  data  set  for  the  purpose  of  comparing  with  a  version  using  200  mb 
steering  since  200  mb  data  were  available  only  since  1962.  A  run 
similar  to  C5-D-PP(s)/(62-81/BT)  using  200  mb  steering  predictors 
instead  of  500  mb  steering  predictors  is  shown  in  Table  5.3  as  C2-D- 
PP(s)/(62-81/BT) 

This  decision  to  limit  the  development  data  to  those  after  1961 
turned  out  to  be  very  fortuitous  since  the  results  indicate  that  the 
largest  overall  decrease  in  error  occurred  between  model  runs 
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C5-D-PP(s)/(46-78/BT)  and  C5-D-PP(s)/(62-81/BT) .  This  surprising  result 
indicates  that  either  the  best  track  cyclone  positions  or  the  analysis 
fields  or  both  for  the  early  pre-satellite  years  are  not  up  to  the 
quality  of  those  since  the  early  60 's  or  that  the  statistical 
characteristics  of  the  analyses  suffered  a  drastic  change  somewhere 
between  the  early  and  later  years. 

A  likely  source  of  the  difference  in  the  analyses  may  stem  from  the 
implementation  of  the  Cressman  (19S9)  scan  analysis  method  at  the 
National  Meteorological  Center  (NMC)  in  19S9.  Similar  data  base 
inhomogeneities  may  occur  on  analyses  after  1978  due  to  the  Cressman 
analysis  method  being  replaced  by  the  'first-guess'/spectral  analysis  PE 
hemispheric  prognosis  package  in  1979.  Neumann,  si  fil*  (1979)  cites 
substantial  deterioration  in  the  performance  of  NHC  statistical -synoptic 
and  statistical-dynamical  models,  particularly  NHC73,  due  to  changes  in 
the  tropical  analysis  error  characteristics. 

Another  surprising  discovery  occurred  when  the  Independent  set  was 
being  verified.  Because  of  the  small  sample  sizes  in  some  of  the 
categories,  a  version  of  the  model  without  any  speed  stratification  was 
developed  on  the  1962-1981  dependent  set.  This  version,  with  either  200 
mb  or  500  mb  steering  predictors,  appears  as  model  runs  C2-D-PP/(62- 
81/BT)  and  C5-D-PP/(62-81/BT)  respectively  in  Table  5.3.  From  the 
results,  it  is  obvious  that  the  speed  stratification  used  in  this  model 
failed  to  enhance  its  predictive  skill  in  any  of  the  categories.  In 
fact.  Table  5.4,  which  presents  the  results  of  the  verification  on  the 
independent  set,  reveals  the  superior  skill  of  the  model  with  no  speed 
stratification.  Compare,  for  example,  model  runs  C5-I-PP(s)  with 
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C5-I-PP  and  runs  C2-I-PP(s)  with  C2-I-PP.  This  is  especially  true  for 
the  version  with  200  mb  steering  predictors  for  the  north-of-ridge 
category. 

Possible  reasons  for  the  failure  of  the  speed  stratification  were 
examined.  Table  5.5  presents  the  number  of  cases  that  occurred  for  a 
given  height  pattern  and  the  subsequent  24-hour  cyclone  speed.  As 
discussed  in  Chapter  3.  the  height  patterns  for  the  Atlantic  were 
computed  according  to: 

H15.40  +  H15.25  "  2  x  H15.-5  for  north-of-ridge 

H10, 10  ■  H10, -15  •  20  for  on-the-ridge 

2  x  H20,5  "  H20.-25  "  H20,25  for  a°uth-of-ridge 

where  ^  is  the  500  mb  height  in  meters  and  the  subscripts  i  and  j 
indicate  the  degrees  of  latitude  and  longitude  away  from  the  cyclone 
center  where  the  heights  are  to  be  observed.  Latitudes  are  positive 
poleward  and  longitudes  are  positive  eastward  from  the  cyclone  center. 
Recall  that  if  the  height  pattern  indicator  HP  was  greater  than  or  equal 
to  zero,  the  cyclone  was  placed  in  the  fast  category  and  if  HP  was  less 
than  zero,  the  cyclone  was  placed  in  the  slow  category.  Although  the 
stratification  worked  well  for  the  on-the-ridge  category,  it  did  poorly 
for  the  other  two  categories.  For  a  given  value  of  the  height  pattern 
Indicator,  there  are  as  many  cases  that  were  fast  as  there  were  slow. 

The  scheme  did  a  poor  job  of  distinguishing  the  fast  cases  from  the  slow 
cases  for  the  given  height  pattern  indicator.  The  results  in  Table  5.4 
bear  this  out. 

Similar  speed  stratification  problems  occurred  in  the  NW  Pacific  as 
shown  in  Table  5.6.  Although  it  appears  that  the  height  pattern 
Indicator  was  able  to  distinguish  between  the  very  fast  and  the  very 


TABLE  5.5 


Number  of  Atlantic  cases  that  occurred  for  a  given  magnitude  of  500  mb 
height  pattern  indicator  (in  meters)  and  subsequent  24-hour  cyclone 
speed  (in  m/s). 


S  <  2.5 

2.5  i  S  i  7.5 

S  >  7.5 

North  78 

238 

144 

100  <  HP 

IS 

59 

46 

50  <  HP  1  100 

12 

95 

37 

0  i  HP  i  50 

24 

68 

33 

-50  i  HP  <  0 

22 

101 

43 

-100  1  HP  <  -50 

175 

408 

150 

HP  <  -100 

Qn  4 

27 

17 

100  <  HP 

4 

51 

18 

50  <  HP  i  100 

26 

132 

22 

0  i  HP  1  50 

39 
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2 

-50  1  HP  <  0 

15 

48 

1 

-100  1  HP  <  -50 
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14 

0 

HP  <  -100 

South  48 

242 

72 

100  <  HP 

16 

112 

21 
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10 
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0  i  HP!  50 

10 

84 

34 

-50  I  HP  <  0 

16 

140 

14 

-100  i  HP  <  -50 

74 

288 

9 

HP  <  -100 

slow  cases  based  on  strongly  positive  or  strongly  negative  values  of  HP, 
this  distinction  is  not  as  clear  as  the  value  of  HP  approaches  zero. 

For  example,  when  the  height  pattern  indicator  was  greater  than  100 
for  the  north-of-ridge  category,  the  scheme  was  able  to  correctly  select 
almost  all  of  the  extremely  fast  cases.  Only  25  cases  that  actually 
were  slow  cases  were  improperly  selected.  Similarly,  the  scheme 
performed  creditably  for  the  on-the-ridge  category  when  the  height 
pattern  Indicator  was  greater  than  100  and  for  the  south-of-ridge 
category  when  a  negative  height  pattern  indicator  was  used  to  select  the 
extremely  slow  cases. 

The  scheme  did  poorly,  however,  on  all  other  occasions.  The  worst 
examples  are  for  the  on-the-ridge  and  south-of-ridge  categories  when  the 
height  pattern  indicator  was  between  zero  and  100.  More  of  the  cases 
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selected  for  the  fast  set  actually  turned  out  to  be  slow  cases. 

Perhaps,  the  speed  stratification  could  have  been  made  on  very  large 
positive  or  negative  values  of  HP  instead  of  zero  as  was  done  in  the 
stratification  here.  But  this  would  have  resulted  in  a  very  small 
sample  size.  Further,  the  fast  and  slow  equations  could  have  been 
applied  to  only  those  cyclone  cases  with  initial  speeds  greater  than  7.5 
m/s  or  less  than  2.5  m/s,  respectively,  as  was  done  in  a  similar  manner 
for  NHC-67.  Neither  approach  was  pursued  here  and  the  decision  was  made 
to  use  the  regression  equations  without  any  speed  stratification.  All 


further  reference  to  CSU84  for  the  Atlantic  will  refer  to  the  version 


without  any  speed  stratification 


Table  5.7  compares  the  skill  of  CSU84  to  the  official  forecast 
error  and  the  operational  error  of  the  best  objective  technique  for  a 
nearly  homogeneous  sample  for  the  period  1973-1979.  Direct  comparison 
between  the  operational  and  CSU84  errors  are  difficult  since  the  CSU84 
cases  during  the  period  are  a  combination  of  dependent  and  independent 
cases.  For  this  comparison,  just  the  dependent  and  independent  cases 
for  1973-1979  from  the  larger  1962-1981  CSU84  runs  were  used.  The  on- 
the-ridge  and  north-of-ridge  categories  were  combined  to  form  the  north 
zone.  This  appeared  reasonable  since  most  of  the  on-the-ridge  cases 
were  north  of  24.5°N. 

These  results  for  CSU84  are  based  on  perfect-prog  data,  best  track 
initial  positions  and  200  mb  steering  predictors  and  appear  to  indicate 
that  the  accuracy  of  CSU84  is  comparable  to  the  best  operational 
technique  at  the  short  forecast  period  and  much  better  at  the  longer 
forecast  periods.  Relative  to  the  other  objective  techniques  currently 
used  by  the  National  Hurricane  Center,  CSU84  appears  far  superior.  For 
example,  the  72-hour  forecast  errors  for  the  north  zone  range  from  401  n 
mi  for  NHC72  to  499  n  mi  for  NHC67  and  for  the  south  zone,  the  errors 
range  from  309  n  mi  for  HURRAN  to  355  n  mi  for  NHC72. 

For  the  purpose  of  comparison  with  a  no-skill  forecast.  Table  5.8 
presents  24-  and  48-hour  forecast  errors  based  only  on  an  extrapolation 
of  the  cyclone's  past  motion  using  the  CSU84  1973-1979  dependent  and 
independent  sets  combined.  The  Table  reveals  that  the  use  of  past  12- 
hour  motion  for  persistence  is  far  better  than  either  the  past  24-hour 
motion  or  the  average  of  the  past  12-  and  past  24-hour  motion.  But  as 
the  results  indicate,  pure  persistence  is  a  poor  indicator  of  future 
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motion  outside  the  deep  tropics.  Note  that  forecast  errors  at  48  hours 
are  well  above  400  n  mi  for  the  north-  and  on-the-ridge  categories 


compared  to  the  48-hour  errors  of  CSU84  that  are  about  half  as  large. 

In  the  development  of  this  new  tropical  cyclone  prediction  model, 
the  emphasis  was  on  improving  the  track  forecast  skill  at  the  higher 
latitudes  for  the  longer  forecast  periods.  These  results  appear  to 
indicate  that  this  goal  has  been  achieved.  Moreover,  the  remarkable 
accuracy  of  the  model  at  24  hours  was  unexpected.  On  reexamination, 
there  are  several  reasons  for  the  model's  very  good  performance  for  both 
the  short-  and  long-term  forecast.  These  reasons  will  be  amplified 
later  in  this  chapter. 

Operational  Simulation.  The  verification  of  this  new  model  thus 
far  has  used  best  track  data  and  a  perfect-prog  approach.  In  order  to 
simulate  more  realistic,  operational  conditions,  the  sensitivity  of  the 
model  was  evaluated  by  introducing  errors  into  the  initial  best  track 
positions.  For  the  seven  year  period  between  1968  and  1974.  Neumann 
(1975)  cites  an  average  initial  position  error  of  25  n  ml,  with  the 
average  error  decreasing  slightly  to  22  n  mi  for  cyclones  within  500  n 
ml  of  the  United  States  (US).  With  this  as  a  guide  and  using  the  model 
run  with  200  mb  steering  for  the  data  period  1962-81  (C2-D-PP(s)/ (62- 
81/BT))  as  the  control  case,  the  initial  position  was  perturbed  by  .3 
degrees  (18  n  ml)  for  those  cyclones  west  of  the  70th  meridian  and  by  .5 
degrees  (30  n  mi)  for  those  cyclones  east  of  the  70th  meridion. 

Since  the  use  of  best  track  Information  to  compute  the  past  12-hour 
and  24-hour  motion  results  in  perfect  persistence  predictors, 
introducing  'erroneous'  initial  positions  can  be  expected  to  result  in 
larger  forecast  errors.  However,  with  the  emphasis  placed  on  this  model 
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for  more  accurate  long-range  position  forecasts,  it  was  hoped  that  the 
deterioration  of  skill  due  to  the  initial  position  error  would  be 
reduced  at  the  longer  time  period.  The  results  of  the  sensitivity  test 
show  that  this  indeed  is  the  case.  Compare  the  results  shown  in  Table 
5.9  with  those  of  C2-D-PP(s)  in  Table  5.3. 

As  one  might  expect,  displacing  the  initial  position  northward  by 
18  to  30  n  mi  degrades  primarily  the  24-hour  forecast  position  and  the 
degradation  decreases  with  time  until  it  is  almost  negligible  at  72 
hours . 

The  degradation  is  most  noticeable  for  cyclones  south  of  the  ridge 
(40%  Increase  in  error  at  24  hours)  and  is  most  pronounced  for  a 
northward  displacement  of  the  initial  position  compared  to  displacements 
southward,  eastward  or  westward.  This  sensitivity  test  indicates  that 
the  use  of  best  track  positions  rather  than  warning  positions  to  produce 
"perfect"  past-motion  predictors  does  not  result  in  an  unrealistically 
strong  reliance  on  the  persistence  component  of  the  model  beyond  24 
hours,  particularly  for  those  cylones  on  and  north  of  the  subtropical 
ridge. 

One  likely  reason  for  the  relatively  small  increase  in  error  due  to 
initial  position  error,  even  for  the  short-term  forecast,  is  that  the 
model  does  not  use  the  cyclone's  current  motion  as  one  of  its 
predictors.  Statistical  models,  such  as  CLIPER  and  NHC73  that  do  use 
current  motion  will  undoubtedly  suffer  a  large  increase  in  forecast 
errors  at  the  shorter  forecast  periods  whenever  large  initial  position 
errors  are  encountered  operationally.  Even  a  relatively  small  initial 
position  error  could  be  disastrous  if  the  error  projects  a  false  change 
in  direction.  Ose  of  past  12-hour  and  past  24-hour  motion  instead  of 
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current  notions  could  be  expected  to  'soften'  the  impact  of  any  initial 
position  error. 

However,  the  use  of  "real"  prognostic  data  in  lieu  of  analysis 
data,  as  was  done  in  this  perfect-prog  approach,  will  almost  certainly 
degrade  the  results;  just  how  much  is  questionable  and  must  be  tested  in 
an  operational  environment.  With  more  sophisticated  numerical  models 
and  more  accurate  prognostic  fields,  the  degradation  will  hopefully  be 
minimal. 

Again,  because  of  the  new  approach  taken  here  of  segmenting  future 
72-hour  cyclone  motion  into  three  discrete  24-hour  time  steps,  the 
dependence  of  the  scheme  on  any  one  prognostic  field  is  limited  to  24 
hours  of  cyclone  movement.  Thus,  the  extent  of  any  erroneous  track 
forecast  based  on  a  "bad"  prognostic  field  should  be  kept  to  a  minimum 
and  hopefully  corrected  on  a  subsequent  24-hour  forecast.  In  addition, 
because  the  design  of  this  model  permits  forecaster  intervention  in 
providing  the  input  data,  if  the  analysis  or  prognostic  field  appears  to 
be  obviously  in  error,  the  input  data  required  by  the  model  could  be 
easily  modified  prior  to  the  computation  of  the  forecast  position. 

Another  desirable  feature  of  this  model  is  that  it  does  not  use 
prognostic  data  to  formulate  the  24-hour  forecast  position.  Nor  does  it 
use  the  current  motion  of  the  cyclone.  The  synoptic  predictors  required 
for  calculating  the  first  24  hours  of  motion  come  from  the  current  and 
24-hour  old  analyses.  Further,  only  the  24-hour  prog  fields  are  used 
for  computing  the  24-  to  48-hour  motion  and  likewise,  only  the  48-hour 
prog  fields  are  used  to  compute  the  48-  to  72-hour  motion.  This  may 
explain  the  sizeable  degradation  in  skill  experienced  by  traditional 
statistical -dynamical  cyclone  prediction  models  that  apply  24-.  36-,  and 


48-hour  forecast  fields  for  each  and  every  forecast  period  in 
transitioning  from  a  perfect-prog  approach  used  in  development  to  an 
operational  setting  that  uses  read  progs. 

Figure  5.1,  taken  from  Neumann  and  Lawrence  (1973),  is  presented  as 
an  example  of  the  predictors  selected  by  the  NHC73  screening  program  for 
zonal  and  meridional  motion  for  zone  22  which  is  centered  near  the 
eastern  tip  of  Cuba.  Their  HOOIO  refers  to  the  current  1000  mb  height, 
H2405  refers  to  the  24-hour  forecast  height  of  the  500  mb  surface,  H4805 
refers  to  the  48-hour  forecast  height  of  the  500  mb  surface,  and  so  on. 
Note  that  their  short  range  12-  and  24-hour  forecasts  are  heavily 
dependent  on  forecast  fields.  Specifically,  10  of  the  12  predictors  for 
24-hour  zonal  motion  come  from  forecast  fields  while  7  of  the  12 
predictors  for  24-hour  meridional  motion  come  from  forecast  fields.  And 
a  fair  number  of  these  predictors  come  from  the  48-hour  forecast  field. 
Since  forecast  fields  for  the  tropics  are  known  to  be  rather  poor  and 
zone  22  can  be  considered  to  lie  in  the  tropics,  the  use  of  forecast 
fields,  particularly  extended  forecasts  for  36  and  48  hours,  can  be 
expected  to  play  a  significant  role  in  the  degradation  of  statistical 
forecast  models  when  actual  forecast  fields  are  used  operationally. 

Model  Runs  Using  Operational  Data.  To  determine  the  amount  of 
degradation  that  might  occur  in  using  this  newly  developed  model 
operationally,  the  model  was  tested  using  actual  forecast  fields  from  a 
numerical  model.  This  was  accomplished  using  prognostic  fields  archived 
at  the  Fleet  Numerical  Oceanography  Center  (FN0C)  for  the  period  1968- 
1981.  Twenty-four-  and  48-hour  forecast  fields  of  sea  level  pressures 
and  500  mb  D-values  on  a  63x63  hemispheric  grid  were  used  to  Introduce 


Fig.  5.1.  Predictors  selected  by  the  NHC73  screening  program  for  zonal  1 

and  meridional  motion  for  zone  22.  H0005  refers  to  the  | 

current  500  mb  height,  H3605  refers  to  the  36-hour  forecast  1 

height  of  the  500  mb  surface,  H0007  refers  to  the  current  j 

700  mb  height,  H0010  refers  to  the  current  1000  mb  height  and  1 

so  on.  The  tables  are  taken  from  Neumann  and  Lawrence  (1973).  i 
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synoptic  data  at  the  appropriate  time.  A  standard  height  of  SS74  m  was 
used  to  convert  the  500  mb  D-values  to  geopotential  heights. 

Unfortunately.  200  mb  prognostic  height  fields  are  not  archived  by 
the  FNOC.  Comparative  testing  is  thus  limited  to  the  version  of  the 
model  using  500  mb  steering  predictors  instead  of  200  mb  steering 
predictors.  Although  comparing  the  independent  set  results  of  model  run 
C5-I-PP/ ( 62-81 /BT) )  and  C2-I-PP/(62-81/BT))  in  Table  5.4  would  indicate 
that  the  version  with  500  mb  steering  is  not  as  accurate  as  the  version 
with  200  mb  steering,  the  verification  runs  using  actual  forecast  fields 
should  be  able  to  Indicate  the  amount  of  degradation  in  forecast 
accuracy  that  the  model  would  suffer. 

Table  5.10A  presents  a  comparison  of  the  model  using  perfect-prog 
data  (C5-DI-PP/ (68-81/BT) )  and  actual  forecast  fields  (C5-DI-AP/ (68- 
81/BT) ) .  Verification  of  the  model  using  both  actual  forecast  fields 
and  perturbation  of  the  initial  position  (C5-DI-AP/ (68-81/PBT) )  is  also 
shown.  To  simulate  warning  positions,  the  best  track  initial  positions 
were  randomly  perturbed  by  an  amount  equal  to  warning  position  er> 
observed  operationally  (Neumann.  1975).  The  relatively  small  amount  of 
degradation  resulting  from  the  use  of  actual  forecast  fields  and 
pseudo-warning  positions  is  comforting.  As  expected,  the  largest 
increases  in  error  using  real  progs  occur  for  the  south-of-ridge  and 
on-the  ridge  categories.  The  poor  quality  of  forecast  fields  in  the 
tropics  is  well-known.  This  is  precisely  where  statistical  schemes  that 
do  not  rely  on  synoptic  data  perform  well.  Yet,  the  results  indicate 
that  the  skill  of  this  model  compares  very  favorably  with  these  analog 
techniques  that  use  just  climatology  and  persistence  and,  as  will  be 


shown  later,  outperforms  them  at  the  longer  forecast  periods  for  the 
south  zone.  At  higher  latitudes,  this  model  appears  superior. 

The  forecast  errors  for  the  on-the-ridge  category  are  the  key  to 
whether  the  model  is  capable  of  discerning  a  major  track  change  from  one 
with  a  westward  component  to  one  with  an  eastward  component.  Although 
there  appears  to  be  a  sizeable  degradation  in  skill  by  introducing 
actual  forecast  fields,  bear  in  mind  that  model  run  C5-DI-PP  was  an 
attempt  to  compare  homogeneous  cyclone  cases  for  the  period  1968-1981. 

By  keeping  the  sample  cases  homogeneous,  however,  both  dependent 
and  independent  cases  from  the  perfect-prog  approach  were  combined. 

Table  5.10B  shows  the  results  from  the  model  run  (C5-I-PP)  on  just  the 
independent  set  using  perfect-progs  for  the  period  1968-1981.  Using 
these  forecast  errors  for  oomparlson,  the  degradation  from  using  actual 
prognostic  fields  is  relatively  small.  A  truly  homogeneous  sample  would 
probably  have  resulted  in  forecast  errors  somewhere  between  the  two 
model  runs  of  C5-DI-PP  shown  in  Table  5.10A  and  model  run  C5-I-PP  shown 
in  Table  5.10B.  Be  aware  that  although  the  1968-1981  data  set  is  a 
mixture  of  dependent  and  Independent  data  for  the  perfect-prog  case,  it 
can  be  considered  totally  Independent  beyond  24  hours  for  the  actual 
prog  case.  From  these  results,  it  would  appear  safe  to  assume  that  the 
use  of  this  model  in  an  operational  setting  would  not  degrade  its 
forecast  skill  as  much  as  the  deterioration  suffered  by  the  other 
statistical -dynamical  models.  Possible  explanations  why  the  approach 
used  in  the  development  of  CSU84  limits  the  normal  operational 
degradation  suffered  by  statistical  models  will  be  offered  in  the 
summary  to  this  chapter. 


5.2  NW  Pacific  and  South  China  Sea 

The  development  of  regression  equations  for  the  NW  Pacific  and 
South  China  Sea  followed  along  the  same  line  taken  for  the  Atlantic.  As 
was  shown  in  Table  5.6,  speed  stratification  was  not  beneficial  for  the 
Pacific;  thus,  only  the  versions  of  the  model  with  no  speed 
stratification  will  be  discussed  from  here  on.  Table  5.11  summarizes 
the  various  model  runs  for  the  NW  Pacific  for  easy  reference.  Table 
5.12  presents  verification  results  for  the  dependent  and  independent 
sets  using  either  the  500  mb  steering  predictors  or  200  mb  steering 
predictors.  Unlike  the  Atlantic  where  the  version  of  the  model  using 
200  mb  steering  was  more  accurate  overall,  the  version  with  500  mb 
steering  was  more  accurate  overall  in  the  NW  Pacific,  except  for  intense 
cyclones  which  are  accelerating  north  of  the  subtropical  ridge.  The 
reason  for  this  is  unclear  although  a  likely  explanation  is  the 
difference  in  the  percent  of  cases  in  the  north-of-ridge  category 
compared  to  the  cases  in  the  south-of-ridge  category  for  the  two  ocean 
basins.  As  will  be  shown  later,  even  for  the  NW  Pacific,  the  version 
using  200  mb  steering  was  more  accurate  for  those  strong  cyclone  cases 
undergoing  recurvature  north  of  the  ridge.  But  since  the  number  of 
these  cases  is  small  compared  to  the  other  cases  for  the  Pacific,  the 
accuracy  of  the  version  with  500  mb  steering  appears  better  overall. 

The  model's  better  performance  on  strong  cyclones  undergoing 
recurvature  when  200  mb  steering  is  used  instead  of  500  mb  steering 
appears  to  support  recent  findings  correlating  the  recurvature 
phenomenon  with  upper  tropospheric  winds.  A  number  of  these  studies 
were  received  earlier  in  Chapter  2. 
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TABLE  5.11  (cont’d) 


In  the  Atlantic,  however,  the  percent  of  cases  in  the  north-of- 
ridge  category  is  rather  high  compared  to  the  number  of  cases  in  the 
other  categories.  Thus,  the  skill  of  the  model  with  200  mb  steering  for 
the  north-of-ridge  category  translates  to  better  skill  overall. 

Table  5.12  also  shows  the  verification  results  for  the  South  China 
Sea.  The  verification  for  the  dependent  set  uses  the  regression 
equations  that  were  developed  using  the  surface  pressure,  850  mb 
heights,  700  mb  heights,  and  500  mb  heights  with  no  ridge  or  speed 
stratification.  The  independent  set  is  comprised  of  the  same  cyclone 
cases  but  uses  the  regression  equations  developed  for  the  NW  Pacific 
cyclones.  Note  that  comparable  results  are  obtained  using  the  NW 
Pacific  regression  equations.  Thus,  unless  it  is  impossible  to  select 
the  proper  ridge  category,  it  is  suggested  that  the  regression  equations 
developed  for  the  NW  Pacific  be  used  for  the  South  China  Sea. 

Table  5.13  presents  a  comparison  between  the  CSU84  version  with  500 
mb  steering  and  perfect-prog  data  and  the  official  forecast  and  the  best 
objective  technique  for  the  period  1979-1981.  These  operational  errors 
were  compiled  from  the  Annual  Tropical  Cyclone  Report  (1979-1981).  Over 
the  last  few  years,  the  best  objective  technique  has  been  the  One-Way 
Interactive  Tropical  Cyclone  Model  (0TCM) .  Although  one  of  the 
objectives  in  the  development  of  a  new  prediction  model  was  the 
Improvement  at  the  longer  forecast  periods,  sizeable  reduction  in 
forecast  error  is  also  evident  at  the  shorter  period. 

Analogous  to  the  Atlantic,  the  sensitivity  of  the  model  was 
evaluated  by  perturbing  the  initial  cyclone  position  as  a  way  of 
simulating  operational  conditions.  Using  the  mean  deviation  of 
satellite-derived  cyclone  positions  from  the  JTWC  best  track  positions 


TABLE  5.13 
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for  1974-1983  (calculated  from  figures  in  Table  2-3  of  the  1983  Annual 
Tropical  Cyclone  Report  of  13.3  n  mi  for  cyclones  possessing  a  well- 
defined  eye,  19.3  n  mi  for  cyclones  possessing  a  well-defined 
circulation  center  and  35.8  n  mi  for  all  others),  the  following 
adjustments  to  the  Initial  best  track  positions  were  made: 

1)  if  the  intensity  of  the  cyclone  was  greater  than  69  kts,  the 
initial  position  was  altered  by  .2  degrees  (12  n  mi); 

2)  if  the  intensity  of  the  cyclone  was  between  35  and  69  kts,  the 
initial  position  was  altered  by  .3  degrees  (18  n  mi); 

3)  otherwise,  the  Initial  position  was  altered  by  .5  degrees 
(30  n  mi); 

Based  on  the  sensitivity  test  done  on  the  Atlantic  cases  indicating 
that  initial  position  errors  increase  the  forecast  errors  the  most  when 
the  estimated  Initial  position  is  in  error  to  the  north  of  the  actual 
position,  all  of  the  best  track  initial  latitudes  were  Increased 
northward  by  the  adjustments  cited  above.  Model  run  C5-I-PP/(62-81/BT) 
using  500  mb  steering  on  the  1962-1981  independent  set  was  used  as  the 
control  run.  The  results  shown  in  Table  5.14  are,  therefore,  for  the 
worst  case  situation.  In  addition,  since  cyclone  intensities  were  not 
recorded  on  every  forecast  situation  prior  to  1977,  and  since  these 
cases  appear  with  an  intensity  of  zero  in  the  best  track  file,  the 
initial  positions  were  altered  by  an  amount  greater  than  necessary  on 
many  of  the  cases.  Also,  with  the  1983  typhoon  season,  JTWC  no  longer 
extrapolates  forward  2-3  hours  to  estimate  a  warning  position.  Instead, 
the  reconnaissance  fix  is  made  at  the  synoptic  time  of  the  warning  and 
no  extrapolation  of  the  initial  warning  position  is  necessary.  Thus, 
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TABLE  5.14 

Forecast  errors  (n  ml)  In  the  NW  Pacific  and  the  percent  increase  in 
error  as  a  result  of  displacing  the  best  track  initial  position 
northward  (PBT)  either  12,  18  or  30  n  ml  depending  on  cyclone  Intensity. 
Model  run  C5-I-PP/ (62-81/BT)  was  used  as  the  oontrol  case  and  appears  as 
BT  in  the  table. 


24-HR 

48-HR 

72-HR 

£1 

PBT 

& 

El 

PBT 

& 

El 

PBT 

S» 

North 

95(75) 

113(75) 

19% 

171(57) 

192(57) 

12% 

242(42) 

269(42) 

11% 

fin 

81(131) 

97(131) 

20% 

161(105) 

185(105) 

15% 

252(91) 

279(91) 

11% 

80(667) 

95(667) 

19% 

155(617) 

167(617) 

8% 

215(570) 

224(570) 

4% 

All 

81(873) 

97(873) 

157(779) 

171(779) 

9% 

221(703) 

234(703) 

6% 

initial  position  error  now  can  be  expected  to  be  reduced  somewhat  from 
those  prior  to  1983. 

As  with  the  Atlantic  sensitivity  test.  Table  5.14  indicates  that 
the  impact  of  a  poor  initial  position  is  greatest  at  24  hours  and 
decreases  steadily  to  about  an  11%  Increase  in  error  at  the  higher 
latitude  at  72  hours.  Unlike  the  Atlantic,  however,  the  effect  of 
displacing  the  initial  position  northward  was  smallest  south  of  the 
ridge. 

Model  Runs  Using  Operational  Data.  The  model  was  tested  on  NW 
Pacific  cyclones  using  actual  forecast  fields  produced  by  the  fleet 
Numerical  Oceanography  £enter  (FNOC)  hemispheric  primitive  equation 
model.  Again,  because  200  mb  prognostic  heights  are  not  archived,  the 
version  of  the  model  with  500  mb  steering  was  used.  For  the  NW  Pacific, 
this  presented  no  problems  since  this  version  was  more  accurate  overall 
than  the  version  of  the  model  with  200  mb  steering.  Table  5.15A 


presents  the  verfication  of  the  model  using  actual  forecast  fields  and 
pseudo-warning  positions  (C5-DI-AP/ (68-81/PBT) )  along  with  the  perfect- 
prog  version  (C5-DI-PP/(68-81/BT) )  and  the  atual  prog  version  using  best 
track  initial  positions  (C5-DI-AP/ (68-81/BT) ) .  To  simulate  warning 
positions,  pseudo-warning  positions  were  computed  by  randomly  perturbing 
the  best  track  Initial  cyclone  positions  by  an  amount  equal  to  warning 
position  errors  encountered  operationally.  Thus,  if  the  cyclone 
intensity  was  greater  than  69  kts,  the  initial  position  was  perturbed  by 
0.2  degrees  (12  n  mi),  for  intensities  of  35  kts  to  69  kts,  the  initial 
position  was  perturbed  by  0.3  degrees  (18  n  mi),  for  intensities  less 
than  35  kts,  the  initial  position  was  perturbed  by  0.5  degrees  (30  n 
mi). 

As  in  the  Atlantic  simulation,  model  run  C5-DI-PP  was  an  attempt  to 
compare  homogeneous  cyclone  cases  for  the  period  1968-1981  and  as  a 
result,  both  dependent  and  Independent  cases  from  the  perfect-prog 
approach  were  combined.  Table  5.15B  presents  the  forecast  errors  using 
just  the  Independent  cases  for  this  time  period.  A  fair  comparison 
between  the  perfect-prog  and  real  prog  verification  should  lie  somewhere 
between  the  two  runs.  Based  on  this  comparison,  the  increase  in 
forecast  error  using  actual  forecast  fields  ranges  from  approximately 
37%  at  72  hours  for  the  north-of-ridge  category  to  approximately  13%  at 
72  hours  for  the  on- th e-ridge  category. 

A  homogeneous  comparison  of  selected  techniques  for  1979-1981 
accomplished  by  T.  Tsui  (personal  communication)  is  shown  in  Table  5.16. 
For  comparison,  the  forecast  errors  for  CSU84  for  the  same  years  with 
actual  forecast  fields  and  with  the  best  track  initial  positions 
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TABLE  5.16 

A  nearly  homogeneous  comparison  of  forecast  errors  (n  mi)  between  CSU84 
using  500  mb  steering  with  actual  forecast  fields  and  perturbed  initial 
position  (C5-DI-AP)  and  selected  dynamical  and  climatological 
operational  objective  techniques  for  1979-81  in  the  NW  Pacific  (provided 
by  Tsui,  NEPRF).  In  each  case,  the  corresponding  official  forecast 
error  is  indicated  in  parentheses  for  a  homogeneous  sample.  HPAC  -  half 
persistence  and  climatology;  CLMW  -  climatology  based  on  warming 
position;  OTCM  -  one-way  interactive  dynamical  model;  NTCM  -  nested 
two-way  interactive  dynamical  model. 


m 

72h 

OTCM 

126(124) 

229(233) 

320(348) 

HPAC 

137(126) 

231(229) 

320(339) 

CLMW 

157(126) 

247(229) 

327(343) 

NTCM 

143(124) 

248(235) 

354(368) 

C5-DI-AP 

80(129) 

152(235) 

235(330) 

randomly  perturbed  (C5-DI-AP/(79-81/PBT) )  are  also  included.  Note  the 
very  substantial  reduction  in  the  forecast  error. 

Since  the  largest  forecast  errors  are  normally  associated  with 
recurving  cyclones,  an  in-depth  analysis  of  the  forecast  errors  for  all 
recurving  cyclones  during  the  period  1978-1981  was  accomplished.  Table 
5.17  and  Table  5.18  list  all  of  the  cyclones  that  comprise  this  set 
along  with  their  respective  official  forecast  errors.  Note  that  with 
the  exception  of  1981,  the  average  error  for  recurving  cyclones  was 
significantly  greater  at  72  hours  than  the  average  error  for  non¬ 
recurving  cyclones. 

The  skill  of  any  objective  technique  and  the  confidence  a 
forecaster  places  on  a  particular  technique  is  based  on  its  ability  to 
forecast  a  change  in  the  track,  either  in  direction  or  epeed.  As 
mentioned  earlier,  any  simple  technique  based  on  persistence  and 
climatology  performs  satisfactorily  on  straight-line  climatological 
tracks.  Unfortunately,  these  objective  aids  do  poorly  on  cyclones  that 


deviate  from  their  previous  course  or  from  climatology.  Thus,  most 
objective  techniques  are  least  useful  in  those  particular  situations 
when  the  forecaster  needs  objective  guidance  the  most. 

To  evaluate  how  the  new  model  performed  on  those  difficult  cases 
when  the  forecaster  really  needed  to  rely  on  objective  guidance,  the  00Z 
and  12Z  forecast  situations  that  resulted  in  abnormally  high  forecast 
errors  (greater  than  500  n  mi.  at  72  hours)  were  selected  from  the 
1978-1981  set  of  recurving  cyclones  presented  in  Tables  5.17  and  5.18. 
These  forecast  situations,  along  with  their  subsequent  forecast 
verification  are  presented  in  Tables  5.19  through  5.22.  The  tables  also 
inlude  a  homogeneous  comparison  of  all  forecasts  made  for  each  cyclone. 
In  comparing  the  official  forecast  errors  and  this  model's  forecast 
errors,  several  points  should  be  kept  in  mind.  The  official  forecast  is 
highly  subjective  and  is  strongly  affected  by  political  considerations 
and  heavily  biased  towards  the  previous  forecasts.  In  contrast,  all 
prediction  models,  including  this  model,  are  objective  techniques  with 
little  continuity  from  one  forecast  to  the  next.  In  these  respects,  the 
official  forecast  tends  to  lag  the  objective  techniques  when  a  large 
change  in  direction  becomes  necessary  and  thus  is  at  a  disadvantage 
whenever  forecast  errors  are  compared.  In  addition,  prior  to  the  1983 
tropical  cyclone  season,  the  initial  warning  position  was  extrapolated 
forward  approximately  three  hours  from  the  latest  reconnaissance  fix 


TABLE  5.17 


List  of  official  forecast  errors  for  NW  Pacific  recurving  cyclones  for 


1978  and  1979.  Number 

of  forecasts  are 

in  parentheses. 

TY  Indicates 

typhoon;  TS  indicates 

tropical  storm. 

im 

2i=flB.-By,Pr 

48-HR  Error 

7  2- HR  Error 

#1  TS  Nadine 

185(16) 

568(12) 

980(08) 

#2  TY  Olive 

100(32) 

224(28) 

328(20) 

#3  TS  Polly 

93(12) 

139(08) 

208(02) 

#7  TY  Virginia 

112(39) 

231(35) 

399(31) 

#15  TY  Faye 

158(40) 

360(36) 

514(29) 

#18  TY  Irma 

92(12) 

134(08) 

154(04) 

#19  TY  Judy 

127(14) 

242(10) 

346(06) 

#22  TY  Mamie 

182(14) 

386(10) 

722(06) 

#24  TY  Ora 

124(17) 

314(11) 

460(05) 

#27  TY  Phyllis 

132(24) 

263(20) 

436(15) 

#29  TS  Tess 

108(16) 

194(12) 

367(08) 

#31  TY  Viola 

96(25) 

269(21) 

434(17) 

#32  TS  Winnie 

238(08) 

614(04) 

— 

Avg.  for  recurvers 

129(269) 

287(215) 

451(151) 

Avg.  for  non-recurvers 

125(305) 

255(220) 

370(153) 

1979 

#2  TY  Bess 

114(17) 

265(13) 

348(09) 

#3  TY  Cecil 

87(37) 

191(33) 

320(29) 

#4  TS  Dot 

130(23) 

244(20) 

315(16) 

#12  TY  Irving 

163(34) 

286(30) 

441(28) 

#13  TY  Judy 

105(36) 

173(27) 

277(23) 

#15  TS  Ken 

116(10) 

278(07) 

415(03) 

#16  TY  Lola 

88(21) 

172(19) 

287(14) 

#19  TY  Owen 

146(33) 

250(29) 

327(25) 

#21  TS  Roger 

195(13) 

251(09) 

303(04) 

#23  TY  Tip 

135(56) 

259(52) 

345(48) 

#27  TY  Abby 

AM(48) 

286(39) 

2311161 

Avg.  for  recurvers 
Avg.  for  non-recurvers 


132(328) 

114(263) 


242(278) 

203(193) 


338(223) 

282(145) 


TABLE  5.18 


List  of  official  forecast 

errors  for  NW 

Pacific  recurving  cyclones  for 

1980  and  1981.  Number  of  forecasts  are 
typhoon;  TS  indicates  tropical  storm. 

in  parentheses. 

TY  Indicates 

1980 

24-HR  Error 

48-HR  Error 

72-HR  Error 

# 2  TS  Carmen 

154(09) 

266(07) 

250(05) 

#3  TY  Dom 

137(39) 

191(27) 

324(23) 

#4  TY  Ellen 

130(31) 

300(27) 

484(23) 

#12  TY  Lex 

137(32) 

314(24) 

499(20) 

#13  TY  Marge 

114(26) 

276(20) 

506(12) 

#17  TY  Orchid 

95(16) 

175(12) 

284(08) 

#20  TY  Sperry 

176(19) 

324(10) 

571(08) 

#21  TS  Thelma 

145(11) 

358(07) 

978(04) 

#22  TY  Vernon 

145(21) 

216(17) 

248(13) 

#23  TY  Wynne 

119(41) 

248(36) 

370(30) 

#25  TY  Betty 

131(36) 

306(29) 

524(28) 

#27  TY  Dinah 

145  (13J 

304(09) 

smon. 

Avg.  for  recurvers 

133(294) 

268(225) 

442(179) 

Avg.  for  non-recurvers 

116(188) 

205(139) 

283(87) 

1981 

#1  TY  Freda 

106(19) 

222(14) 

369(09) 

#2  TS  Gerald 

161(15) 

289(11) 

426(07) 

#5  TY  June 

119(18) 

227(13) 

196(05) 

#18  TY  Agnes 

104(21) 

167(17) 

244(12) 

#21  TY  Doyle 

149(10) 

269(06) 

494(02) 

#22  TY  Elsie 

97(27) 

213(23) 

377(19) 

#24  TY  Gay 

163(32) 

275(27) 

410(24) 

#25  TY  Hazen 

130(33) 

263(30) 

361(26) 

#26  TY  Irma 

1M291 

H8I21L 

1411211 

Avg.  for  recurvers 

120(204) 

222(166) 

324(125) 

Avg.  for  non-recurvers 

122(254) 

216(177) 

342(122) 
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position.  As  a  result,  the  ensuing  initial  position  error  could  be 
expected  to  have  an  impact  on  the  subsequent  forecast  positions. 

However,  as  discussed  earlier,  sensitivity  tests  on  the  model 
reveal  that  initial  position  errors  affect  the  short-term  forecast 
positions  as  one  might  expect  and  the  effect  becomes  almost  negligible 
at  the  longer  forecast  periods. 

Note  that  with  the  exception  of  TS  Roger  in  1979  and  one  forecast 
for  TY  Betty  in  1980,  both  versions  of  CSU84  using  perfect-progs  and  500 
mb  steering/200  mb  steering  consistently  produced  forecast  errors  that 
were  much  lower  than  those  of  the  official  forecasts.  Although  the 
forecast  errors  for  C5-I-AP  (that  simulates  operational  conditions)  were 
usually  greater  than  the  perfect-prog  versions,  the  errors  were 
consistently  less  than  the  official  forecast  errors.  These  results 
appear  to  Indicate  that  this  new  prediction  model  is  able  to  accurately 
predict  cyclone  motion  even  during  very  difficult  situations  that 
resulted  in  official  forecast  "busts". 

The  only  Instances  where  the  model  did  poorly  (but  still  much 
better  than  the  official  forecasts)  were  on  TY  Faye  in  1978  and  on  the 
last  few  forecasts  for  TY  Gay  in  1981.  Based  on  the  total  forecast 
errors  for  the  period  1978-1981,  CSU84  -  using  actual  forecast  fields 
and  pseudo-warning  positions  -  was  able  to  decrease  the  forecast  error 
of  this  homogeneous  sample  by  more  than  30%  for  all  three  forecast 
periods. 

5.3  North  Indian  Ocean 

The  verification  results  for  the  North  Indian  Ocean  with  the 
perfect-prog  method  are  presented  in  Table  5.23.  Although 
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the  comparison  between  the  model  and  the  operational  forecast  is  on  an 
Inhomogeneous  sample,  the  amount  of  improvement  is  very  substantial. 

The  model  was  not  run  for  this  ocean  basin  using  real  progs.  However, 
assuming  the  decrease  in  accuracy  using  actual  forecast  fields  would  be 
on  the  same  order  as  the  Atlantic  and  the  NW  Pacific  —  20%  at  48  hours 
and  30%  at  72  hours  -  errors  of  176  n  mi  at  48  hours  and  282  n  ml  at  72 
hours  can  be  projected  for  all  cases.  Compared  to  the  best  operational 
techniques,  these  forecast  errors  translate  to  an  improvement  of  32%  at 
48  hours  and  50%  at  72  hours!  At  24  hours,  there  is  no  degradation  due 
to  the  use  of  actual  forecast  fields  since  they  do  not  enter  into  the 
computation  of  the  24-hour  forecast  position.  However,  as  shown  for  the 
Atlantic  and  NW  Pacific,  perturbing  the  best-track  initial  position  to 
simulate  warning  position  errors  does  have  a  strong  impact  on  the  24- 
bour  forecast  position.  Even  using  the  worst  case  scenario  of  a 
northward  track  perturbation  for  all  cases  (which  results  in  a  25% 
decrease  in  accuracy),  the  resulting  improvement  in  the  forecast  error 
at  24  hours  oompared  to  HPAC  is  a  very  creditable  27%. 

5.4  Synopsis  and  Summary 

As  mentioned  earlier,  the  universal  problem  of  all  statistical- 
dynamical  prediction  models  that  use  a  perfect-prog  approach  to 
Introduce  prognostic  data  in  the  development  of  regression  equations  is 
the  degradation  in  skill  they  suffer  in  operational  use  when  actual 
forecast  fields  are  used.  In  fact,  other  than  NHC73,  they  do  not 
attempt  to  forecast  beyond  48  hours.  Table  5.24  illustrates  the  amount 
of  degradation  that  statistical-dynamical  models  suffer  when  they  use 


actual  forecast  fields 


131 


TABLE  S.24 

A  comparison  of  development  and  operational  errors  for  selected 
statistical  techniques.  The  developmental  errors  for  NHC73.  NHC72,  and 
CLIPER  are  from  Neumann  and  Lawrence  (1973).  The  operational  errors  for 
NHC73.  NHC72  and  CLIPER  are  the  entire  sample  errors  from  Neumann  and 
Pelissier  (1981).  The  errors  for  SNT  are  from  Nomoto  £t  si>  (1976)  and 
the  errors  for  Keenan  are  from  Keenan  (1984). 


2± 

-M 

48-HR 

1 

72-HR 

DEV 

OP 

(OP-DEV) 

DEV 

OP  (OP-DEV) 

DEV 

OP  (OP-DEV) 

NHC73 

42 

113 

71 

78 

228 

150 

142 

344 

202 

NHC72 

68 

113 

45 

175 

257 

82 

308 

375 

67 

CSU84A*** 

93 

99 

6 

171 

184 

13 

248 

281 

33 

CSU84P**** 

77 

84 

7 

146 

164 

18 

205 

257 

52 

SNT* 

141 

149 

8 

234 

388 

154 

— 

— 

— 

Keenan** 

83 

99 

16 

135 

222 

87 

— 

— 

— 

•  SNT  is  statistical -dynamical  model  developed  by  Nomoto  £l.  (1976) 
for  the  NW  Pacific  in  the  vicinity  of  Japan. 


••  Keenan  is  statistical-dynamical  model  developed  by  Keenan  (1984)  for 
the  Australian  region. 

•••CSU84A  is  the  Atlantic  version  of  the  new  model  using  S00  mb  steering. 
The  errors  were  on  a  homogeneous  sample  for  the  period  1968-1981  and  taken 
from  Table  3.10A. 

••••CSU84P  is  the  NW  Pacific  version  of  the  new  model  using  500  mb 
steering.  The  errors  were  on  a  homogeneous  sample  for  the  period 
1968-1981  and  taken  from  Table  5. ISA. _ 


Note  the  relatively  small  amount  of  deterioration  associated  with 
CSU84.  Although  the  degradation  suffered  by  NHC72  is  also  comparatively 
small.  NHC72  does  not  use  prog  data. 

This  large  degradation  in  skill  associated  with  statistical- 
dynamical  schemes  is  in  large  part  due  to  the  over-reliance  or 
overweighing  of  the  "perfect"  prognostic  data  at  the  expense  of  the 
analysis  data.  This  problem  was  vividly  shown  earlier  with  the  NHC73 
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predictor  selection  process  that  chose  predictors  from  36-  and  48-hour 
forecast  fields  for  the  short-range  12-  and  24-hour  cyclone  motion.  As 
a  result,  although  the  model  may  produce  very  impressive  results  in 
development,  when  actual  forecast  fields  that  are  obviously  not  as 
accurate  as  the  verifying  analysis  fields  are  used  operationally,  the 
accuracy  deteriorates  rapidly.  In  certain  instances  where  prognostic 
data  are  notoriously  poor,  as  in  the  ocean  region  to  the  west  of 
Australia,  regression  technique  using  just  the  analysis  data  outperforms 
one  using  prognostic  data  (T.  Keenan,  1984). 

The  approach  developed  here  attempts  to  circumvent  this  problem  by 
discretizing  72  hours  of  cyclone  motion  into  three  shorter  24-hour  time 
steps.  This  permits  several  important  processes  to  occur: 

(1)  the  introduction  of  a  revised  persistence  predictor  after  each 
24  hours  of  motion 

(2)  the  opportunity  to  change  the  category  relative  to  the  ridge 
after  each  24  hours  of  motion 

(3)  the  limiting  of  the  reliance  of  the  model  on  any  one 
prognostic  field  to  24  hours  of  motion 

(4)  the  absence  of  any  degradation  for  the  24-hour  forecast  since 
the  first  24  hours  of  motion  is  based  only  on  analysis  fields 

(3)  the  updating  of  prognostic  data  relative  to  a  new  cyclone 
position  —  not  the  initial  cyclone  position. 

The  advantages  of  this  method  over  the  recent  customary  method  of 
applying  prognostic  data  are  obvious.  If  the  prognostic  fields  are  in 
error,  there  will  be  no  impact  on  the  24-hour  forecast  position.  In 
addition,  since  CS1784  does  not  use  current  motion  as  a  predictor,  an 
initial  position  error  will  have  a  much  smaller  effect  on  the  24-hour 
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forecast  position  than  It  will  with,  for  example.  NHC73  that  does  use 
current  motion  as  a  predictor.  If  the  24-hour  forecast  position  is 
accurate,  the  updated  persistence  predictor  as  well  as  the  ridge 
category  in  which  to  place  the  cyclone  should  be  accurate.  Then,  even 
if  the  24-hour  forecast  field  is  in  error,  the  resulting  48-hour 
forecast  position  certainly  will  be  more  accurate  than  if  the  motion 
over  48  hours  was  calculated  from  the  initial  position.  Also  recall 
that  half  of  the  data  needed  to  compute  the  height  and  pressure 
tendencies  at  the  24-hour  forecast  position  to  determine  the  48-hour 
forecast  position  come  from  the  analysis  field.  Thus,  a  solid  forecast 
track  for  the  early  portion  of  the  forecast  period  has  been  created  on 
which  to  build  the  final  24  hours  of  motion  to  determine  the  72-hour 


forecast  position. 


k.  In  order  to  show  skill  in  forecasting 


cyclone  motion,  a  model's  forecast  accuracy  must  exceed  that  achieved  by 
basic  methods  such  as  climatology  and  persistence.  Neumann  and 
Pellssier  (1981)  have  proposed  using  the  CLIPER  model  as  a  kind  of 
normalizer  or  no-skill  model  to  provide  a  benchmark  on  which  to  judge 
the  skills  of  other  models.  They  compute  the  percentage  improvement  or 
deterioration  of  a  forecast  as  to  how  well  the  CLIPER  forecast  has  done. 
This  percentage  is  given  by: 

%  -  100  <Ec  -  ly  /  Ec. 

where  %  is  the  percent  improvement  or  deterioration  over  CLIPER  forecast 
error,  Ec  is  CLIPER  forecast  error  and  Eg  is  the  model  forecast  error. 
Positive  %  indicates  the  model  performed  better  than  CLIPER,  whereas 
negative  %  Indicates  the  model  performed  poorer  than  CLIPER. 


-  A 


V 


Figure  5.2  shows  the  relative  skill  of  NHC73,  which  was  the  best 
overall  performer  among  the  operational  models  at  the  US  National 
Hurricane  Center,  and  CSU84.  The  north-of-ridge  and  on-th e-ridge 
categories  for  CSU84  were  combined  to  form  the  north  zone  since  a 
cursory  examination  of  those  cases  in  the  on-the-ridge  category  revealed 
that  nearly  all  of  the  cases  were  north  of  24.5°N.  This  version  of  the 
model  used  actual  forecast  fields  with  random  perturbation  of  the 
Initial  position  to  simulate  operational  conditions  as  closely  as 
possible  for  cyclone  cases  during  1973-1979. 

Panel  a  presents  the  "all"  cyclone  category.  It  appears  that 
relative  to  CLIPER,  CSU84  possesses  six  times  ss  much  skill  as  NHC73  at 
24  and  72  hours.  It  is  also  able  to  beat  the  official  forecast  by  a 
wide  margin  at  all  forecast  periods. 

The  most  significant  feature  is  shown  on  panel  b  where  CSU84  is 
clearly  superior  to  climatology  and  persistence  at  the  longer  forecast 
periods  even  in  the  south  zone.  While  NHC73  is  degrading  rapidly 
between  48  and  72  hours,  CSU84  is  improving  in  performance. 

Finally,  panel  c  compares  the  relative  performance  on  cyclones  in 
the  northern  zone.  The  strength  of  any  statistical-dynamical  model 
should  emerge  in  this  category  and  indeed  NHC73  and  CSU84  performed 
better  than  CLIPER  by  the  largest  margin  in  this  zone.  The  relative 
skill  of  CSU84  appears  to  be  nearly  twice  that  of  NHC73  for  this 
oategory. 

Figures  5.3  through  S.5  offered  a  graphic  oomparlson  of  the 
relative  skill  of  CSU84  and  NHC73  for  the  three  categories  in  the 
Atlantic.  As  shown  in  panel  a  of  Fig.  5.2  and  again  in  Fig.  5.3,  CSU84 
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in  the  forecast  error  over  CLIPER  is  shown  with  the  percent 
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is  six  times  as  skillful  as  NHC73  relative  to  CL1PER  at  both  the  24-  and 
72-hour  foreast  periods  for  the  entire  sample.  CSU84's  high  skill  at 
the  24-hour  period  comes  primarily  from  its  strong  performance  on  north 
zone  cyclones  (see  Fig.  5. 5)  while  the  skill  at  72  hours  comes  largely 
from  its  performance  on  south  zone  cyclones  (see  Fig.  5.4). 

The  model's  encouraging  performance  on  north  zone  cyclones  appears 
to  be  relatjd  to  its  ability  to  accurately  depict  and  anticipate  changes 
in  the  synoptic  and  steering  forcings  as  a  result  of  the  ridge 
stratification  and  the  24-hour  time  steps.  The  short  time  steps  permit 
the  category  the  cyclone  is  placed  relative  to  the  ridge  to  change  as  it 
becomes  necessary.  The  model's  ability  to  provide  accurate  position 
forecasts  at  the  longer  periods  for  the  south  zone  is  similarly  related 
to  the  segmented  24-hour  time  steps  which  permits  the  reintroduction  of 
updated  persistence  at  each  step. 

A  similar  comparison  for  the  NW  Pacific  was  done  using  HPAC  as  the 
normalizer  or  no-skill  model  to  provide  the  benchmark  for  comparison. 

The  forecast  error  data  used  for  this  comparison  are  from  NEPRF  (T. 

Tsui,  personal  communication)  and  appear  in  Fig.  5.6  and  5.7. 

The  relative  skill  of  CSU84  compared  to  the  two  best  overall 
performers  in  the  NW  Pacific  (OTCM  and  NTCM)  is  considerable  and  the 
apparent  improvement  in  forecast  error  provided  by  CSU84  is  encouraging. 
The  model  was  at  least  four  times  as  skillful  as  the  best  operational 
technique  (OTCM)  at  all  forecast  periods.  As  in  the  Atlantic,  the  ridge 
stratification  and  the  24-hour  time  steps  are  believed  to  provide  a 
method  for  the  model  to  optimumly  formulate  the  current  synoptic  state 
and  its  changes  that  affect  cyclone  movement  over  a  24-hour  time  span. 


Percent  Difference  In  Mean  Forecast  Error 
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Fig.  5.6.  Performance  of  specified  models  relative  to  performance  of 
HPAC  on  NW  Pacific  cyclones.  Sample  is  nearly  homogeneous 
over  years  1979-1981.  Forecast  errors  for  official,  HPAC, 
NTCM  and  OTCM  are  from  Tsui  (personal  communication). 


141 


s* 

« 

*-»  6 

£  -90 

c  e 

w  2  -60 
CE 

g^S-30 

SI  o 

fel|+30 

<  5; 
o  p 

uj  t 

a:  <3  ___ 

O  _90 

lx.  e 

2  -60 
?g-30 

-iz 


ALL  NW  PACIFIC  CASES  (1979-81) 


-57  nm  (42%) 


£  -40 


|2  -2o(-  _ 


1 1  nn^8%)  -|  |  ny%)  | 

JTWC  OTCM  CSU84 

+6  nm  (-4%) 


-79  nm  (34%) 


24  HR 


48  HR 


-2nm(9%)  -2nm(9%)  g 

JTWC  CSU84 

-17  nm(-7%) 


-85  nm  (27%) 


72  HR 


JTWC  NTCM  0  nm 


J  ^  w '  rv  i  v 

ju  H  c 

+30  •  +19  nm  (-6 */»1r^b 
*T*  +34  nm  (-IIV*) 


OTCM  081)84 


Fig.  5.7.  Relative  skill  of  CSU84  for  the  NW  Pacific  oompared  to  OTCM 

and  NTCM  using  HPAC  as  a  no-skill  model  to  provide  a  benchmark. 
Decrease  or  increase  in  the  forecast  error  over  HPAC  is  shown 
with  the  percent  improvement  in  parentheses.  Errors  for  OTCM, 
NTCM  and  JTWC  are  from  NEPRF  (T.  Tsui,  personal  communication). 
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Figures  5.8  and  5.9  for  the  Atlantic  and  NW  Pacific,  summarize  the 
results  of  the  various  attempts  to  simulate  operational  conditions  in 
comparison  with  operational  forecast  models.  The  perfect-prog  runs  for 
CSU84  are  shown  along  with  the  runs  using  actual  progs  as  well  as  the 
runs  with  both  actual  progs  and  pseudo-warning  positions.  Note  that  the 
CSU84  forecast  errors  from  the  operational  simulation  runs  lie  closer  to 
the  perfect-progs  than  to  those  of  the  other  operational  models.  This 
appears  to  strongly  indicate  that  CSU84  is  capable  of  improving  one-  to 
three-day  track  predictions  in  the  operational  arena. 

Based  on  the  CSU84  forecast  errors  for  the  Atlantic  and  NW  Pacific, 
it  would  appear  that  cyclones  in  the  Atlantic  basin  are  more  difficult 
to  forecast.  A  cursory  examination  of  the  number  of  forecast  situations 
in  each  of  the  three  categories  relative  to  the  ridge  indicates  that  in 
the  NW  Pacific,  the  south-of-ridge  cases  dominate  with  a  sizeable  number 
in  the  on-the-ridge  category.  Relatively  few  cases  occur  north  of  the 
ridge.  For  the  Atlantic,  the  situation  is  reversed.  The  number  of 
cases  north  of  the  ridge  is  almost  as  much  as  the  number  to  the  south 
with  relatively  few  cases  on  the  ridge. 

The  reason  for  the  observed  difference  in  cyclone  tracks  for  the 
two  ocean  basins  is  due  to  differences  in  climatology.  Majority  of 
cyclones  in  the  NW  Pacific  form  at  a  lower  latitude  than  those  in  the 
Atlantic.  Hence,  Pacific  cyclones  have  a  tendency  to  track  poleward  for 
greater  distances  prior  to  recurving  eastward  than  cylones  in  the 
Atlantic.  Following  recurvature,  shearing  or  extratropical  transition 
follows  rather  quickly.  Atlantic  cylones,  on  the  other  hand,  come  under 
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Fig.  5.8.  Mean  Atlantic  forecast  error  (n  mi)  for  the  various  CSU84 
runs  compared  to  CLIPER  and  NHC73  at  24,  48  and  72  hours. 
Errors  for  CLIPER  and  NHC73  are  for  1973-1979  entire 
sample  and  are  taken  from  Neumann  and  Pelissier  (1981). 

One  of  the  CSU84  runs  denoted  by  oooooooooo  is  for  1973-1979. 
The  other  runs  are  for  1968-1981. 
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Fig.  5.9.  Mean  NW  Pacific  forecast  error  (n  mi)  for  the  various  CSU84 
runs  compared  to  HPAC,  OTCM  and  NTCM  at  24,  48  and  72  hours. 
Errors  for  HPAC,  OTCM  and  NTCM  are  for  1979-1981  and  were 
provided  by  NEPRF  (T.  Tsui,  personal  communication) .  One  of 
the  CSU84  runs  denoted  by  oooooooooo  is  for  1979-1981.  The 
other  runs  are  for  1968-1981. 


the  influence  of  mid-latitude  westerlies  before  tracking  very  far 
poleward  and  the  portion  of  their  track  that  would  fall  into  the  on- 
ridge  category  is  relatively  small. 

In  addition,  as  shown  by  Xu  and  Gray  (1982)  the  percentage  of 
cyclones  undergoing  fast  (>  7.5  m/s),  slow  (<  2.5  m/s)  or  looping  motion 
is  much  higher  overall  in  the  Atlantic  than  in  the  NW  Pacific.  Figure 
5.10  is  taken  from  the  work  of  Xu  and  Gray  (1982)  to  illustrate  the 
difference.  Due  to  the  large  number  of  fast-moving  cyclones  at  low 
latitudes  and  the  large  number  of  slow-moving  and  looping  cyclones  at 
higher  latitudes,  along  with  the  relatively  small  number  of  cases,  the 
average  forecast  errors  become  inflated  by  these  difficult  forecast 
situations. 

In  view  of  these  difficult  forecast  cases,  the  quite  skillful 
performance  of  CSU84  in  the  Atlantic  is  particularly  noteworthy.  In  the 
NW  Pacific  and  North  Indian  Oceans  where  a  statistical-dynamical  model 
is  not  available,  the  potential  for  a  large  improvement  in  track 
forecast  accuracy  looks  very  promising. 
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Fig.  5.10.  Percentage  of  storms  in  the  Atlantic  and  western  North  Pacific 
which  are  fast  (>  7.5  m/s),  slow  (<2.5  m/s)  or  looping  vs. 
latitude.  (From  Xu  and  Gray,  1982). 


6.  CONCLUSION 


Although  significant  progress  has  been  made  in  developing 
statistical  and  dynamical  models  for  forecasting  the  tracks  of  tropical 
cyclones,  the  forecast  errors  over  the  past  decade  have  shown  little 
improvement.  The  lack  of  observations  over  the  data-sparse  regions  of 
the  tropical  oceans,  and  the  resulting  difficulty  of  accurately 
specifying  the  Initial  state  of  the  atmosphere  for  both  the  dynamical 
tropical  cyclone  models  and  the  global  PE  models,  have  most  often  been 
cited  for  the  lack  of  improvement.  ''Give  me  more  observations  and 
better  analyses  and  I'll  give  you  better  track  forecasts"  has  been  the 
motto  for  the  last  few  years. 

This  adage  of  requiring  better  analyses  is  certainly  true,  but 
another  important  consideration  to  improving  the  accuracy  of  statistical 
track  prediction  models  is  the  proper  application  of  the  existing  data. 
The  need  for  synoptic  data,  especially  in  a  prognostic  sense,  to  be  able 
to  forecast  cyclone  motion  at  the  longer  forecast  periods  and  to  predict 
anomalous  motion  is  well  known.  However,  the  statistical-dynamical 
models  in  existence  that  Incorporate  the  output  from  a  numerical  model 
as  one  set  of  predictors  apply  the  prognostic  data  always  from  the 
initial  position  of  the  cyclone  to  the  desired  foreoast  projection  of 
12,  24,  36.  48,  or  72  hours.  The  Inherent  weakness  of  using  prognostic 
data  in  such  a  manner  was  shown  by  comparison  with  the  results  from  this 
newly  developed  model  using  an  approach  of  introducing  prognostic 


data  relative  to  each  new  cyclone  position  as  if  they  were  analysis 
data. 

The  suggested  new  nethod  of  applying  the  various  climatological, 
persistence  and  synoptic  predictors  is  patterned  after  the  normal 
procedure  followed  by  operational  forecasters  as  they  subjectively 
formulate  the  forecast  track.  This  process  normally  begins  with  an 
examination  of  the  current  and  24-hour  old  analysis  fields  relative  to 
the  initial  position  of  the  cyclone.  Synthesis  of  the  synoptic 
information  combined  with  persistence  produces  a  short-range  24-hour 
forecast  motion.  The  24-hour  forecast  field  is  only  then  Introduced  to 
provide  an  updated  environmental  Influence  on  the  cyclone  relative  to 
the  new  position  and  the  cyclone  is  projected  forward  another  24  hours. 
The  48-hour  forecast  field  is  then  used  to  provide  the  environmental 
forcings  relative  to  the  cyclone's  projected  48-hour  forecast  position. 
New  persistence  Information  is  also  formulated  relative  the  latest 
cyclone  position  and  the  final  24  hours  of  motion  is  projected. 

Although  the  scheme  proposed  here  was  developed  with  an  emphasis  on 
improving  the  forecast  accuracy  for  cyclones  at  the  higher  latitudes  and 
at  the  longer  forecast  periods,  it  performs  at  the  lower  latitudes  and 
at  the  shorter  forecast  period  with  significant  skill.  As  shown  in 
Chapter  5.  the  degradation  in  skill  resulting  from  the  use  of  actual 
forecast  fields  was  minimal  oompared  to  that  suffered  by  the  other 
discussed  statistical-dynamical  models. 

The  Increase  in  forecast  errors  between  the  perfect-prog  and  the 
actual  prog  cases  for  the  three  other  statistical-dynamical  models  that 
were  examined  ranged  from  87  n  mi  to  154  n  mi  at  48  hours  oompared  to 


less  than  20  n  mi  for  CSU84.  What's  more,  two  of  these  models  do  not 
attempt  to  forecast  beyond  48  hours. 

A  comparison  of  the  model  under  a  simulated  operational  setting  on 
a  nearly  homogeneous  sample  in  comparison  with  the  best  operational 
techniques  in  the  Atlantic  and  NW  Pacific  using  climatology  and 
persistence  as  a  normalizer  to  indicate  the  relative  skill  of  the  models 
as  compared  to  a  no-skill  forecast  was  presented  in  Chapter  5.  The 
comparison  indicated  the  strength  of  the  new  scheme  at  both  the  shorter 
24-hour  foreast  period  and  the  longer  72-hour  forecast  period.  For  the 
entire  Atlantic  sample.  CSU84  is  six  times  as  skillful  as  NHC73  relative 
to  CLIPER  at  both  the  short  and  long  forecast  times.  CSU84's  relatively 
high  skill  at  the  24-hour  period  comes  primarily  from  its  strong 
performance  on  north  zone  cyclones  while  the  skill  at  72  hours  comes 
largely  from  its  performance  on  south  zone  cyclones. 

The  model's  encouraging  performance  on  north  zone  cyclones  appears 
to  be  related  to  its  ability  to  acurrately  depict  and  anticipate  changes 
in  the  synoptic  and  steering  forcings  as  a  result  of  the  ridge 
stratification  and  the  24-hour  time  steps.  Its  ability  to  provide 
accurate  position  forecasts  at  the  longer  periods  for  the  south  zone  is 
similarly  related  to  the  segmented  24-hour  time  steps  which  permits  the 
reintroduction  of  persistence  at  each  step. 

Equally  strong  performance  of  CSU84  in  the  NW  Pacific  was  evident. 
The  model  was  at  least  four  times  as  skillful  as  the  best  operational 
technique  (OTCM)  at  all  forecast  periods.  Again,  the  ridge 
stratification  and  the  24-hour  time  steps  are  believed  to  provide  a 
method  for  the  model  to  formulate  the  current  synoptic  state  and  its 
changes  that  affect  cyclone  movement  over  a  24-hour  time  span. 


The  key  to  this  large  improvement  in  accuracy  lies  in  the  unique 
application  of  the  predictors  normally  used  in  statistical  prediction 
models.  Climatology  is  incorporated  via  a  stratification  scheme  based 
on  the  position  of  the  cyclone  relative  to  the  500  mb  subtropical  ridge. 
Each  cyclone  forecast  situation  is  initially  stratified  by  one  of  three 
categories  relative  to  the  500  mb  ridge  —  south-of-ridge,  on-the-ridge, 
or  north-of-ridge  —  and  the  appropriate  regression  equations  are  used 
to  calculate  the  24-hour  displacement.  The  72-hour  forecast  period  is 
segmented  into  three  24-hour  time  steps  and  the  24-hour  displacements 
are  stepped  off  until  the  desired  multi-day  forecast  projection  is 
reached.  This  permits  the  introduction  of  updated  persistence  and 
synoptic  predictors  relative  to  the  new  cyclone  position  as  well  as  the 
opportunity  to  change  the  stratification  category  as  it  becomes 
necessary  based  either  on  the  past  direction  of  motion  or  on  the 
Judgment  of  the  forecaster.  Two  hundred  millibar  geostrophic  winds  5° 
to  20°  northwest  and  north  of  the  cyclone  are  used  as  key  steering 
predictors  in  the  Atlantic  and  in  recurvature  situations  for  strong 
typhoons  in  the  NW  Pacific.  Five  hundred  millibar  steering  predictors 
proved  to  be  more  accurate  for  all  other  situations  in  the  NW  Pacific, 
however. 

An  additional  advantage  of  a  prediction  model  such  as  this  is  the 
ease  with  which  it  can  be  applied  directly  at  operational  locations  and 
aboard  ships  without  a  link  to  a  large  computer  data  base.  The  initial 
analyses  and  the  24-  and  48-hour  forecast  fields  of  sea  level  pressure 
and  500  mb  and  200  mb  geopotential  heights  are  all  that  is  required. 
Further,  since  the  model  permits  forecaster  intervention  in  computing 
the  forecast  positions,  the  operational  forecaster  could  easily  modify 


the  synoptic  predictors  required  for  input  if  the  analysis  or  prognostic 
field  appears  to  be  in  error.  The  forecaster  could  also  make  his  own 
forecast  fields  based  on  subjective  extrapolation  of  large-scale 
synoptic  patterns  in  the  areas  where  predictors  are  required. 

This  discretized  forecast  method  revealed  the  significant 
contribution  to  the  forecast  of  the  steering  predictors  on  24-hour 
motion  at  the  higher  latitude  in  the  Atlantic.  For  the  NW  Pacific, 
steering  predictors  had  a  strong  impact  only  for  east-west  motion  north 
of  the  ridge  and  only  at  the  200  mb  level.  The  motion  of  south-of-ridge 
cyclones  was  dominated  by  persistence  whereas  the  motion  of  on-the-ridge 
category  cyclones  was  strongly  influenced  by  the  synoptic  predictors, 
primarily  the  current  S00  mb  heights.  For  the  North  Indian  Ocean, 
cyclone  motion  south  of  the  ridge  was  dominated  by  persistence. 

However,  cyclone  motion  on  or  north  of  the  ridge  was  strongly  dependent 
on  synoptic  predictors,  primarily  the  24-hour  tendency  of  surface 
pressure. 

Incremental  improvement  in  the  accuracy  of  this  model  can  likely  be 
made  by  several  smaller  modifications.  It  is  felt,  however,  that  the 
major  improvement  was  gained  by  the  special  methodological  approach 
developed  in  this  paper. 

One  obvious  source  of  improvement  is  to  redevelop  the  regression 
equations  for  use  in  the  NW  Pacific  and  North  Indian  Oceans  using 
synoptic  data  from  the  Navy  Fleet  Numerical  Oceanography  Center  (FNOC) 
rather  than  the  NMC  data  used  to  develop  the  CSU-84  equations.  Harrison 
(1982),  while  comparing  the  NTCM  and  MFM,  discovered  considerably 
degraded  performance  for  the  NTCM  when  it  was  initialized  with  NMC 
analysis  rather  than  FNOC  tropical  analysis.  An  initialization 


procedure  that  resulted  in  much  smaller  amplitude  of  the  zonal  wind 
features  in  the  NMC  data  (greatest  in  the  lower  1000-850  mb  levels)  was 
suspected  as  the  primary  reason.  Redevelopment  of  the  regression 
equations  usng  the  FNOC  data  was  the  projected  plan  if  the  degradation 
using  the  Navy  prognostic  fields  had  been  severe.  Since  it  was  not,  the 
plan  to  redevelop  the  equations  was  abandoned.  However,  since  the 
statistical  characteristics  of  the  two  data  sources  differ,  one  would 
expect  better  model  performance  if  it  used  the  same  data  source  in 
operational  use  that  it  used  in  development.  The  fact  that  the 
deterioration  in  the  accuracy  of  the  model  was  as  small  as  it  was  is 
certainly  encouraging. 

The  ideal  method  of  developing  the  regression  equations,  as  alluded 
to  in  Chapter  2,  is  to  use  the  MOS  approach  and  use  the  actual  forecast 
fields  instead  of  analysis  fields.  This  would  certainly  be  possible  for 
the  NW  Pacific  since  there  are  adequate  cyclone  cases  and  the  forecast 
fields  for  the  surface  and  500  mb  level  are  available  since  1968.  The 
degradation  in  skill  from  the  development  to  the  operational  fields  then 
primarily  would  be  eliminated.  This  may  now  be  more  difficult,  however, 
since  the  Navy  has  abandoned  its  hemispheric  PE  model  and  is  now  using  a 
global  spectral  model. 

A  second  source  of  possible  improvement  to  this  scheme  would  be  the 
substratification  of  the  development  set  to  separate  those  cyclones  with 
very  fast  motion  07.5  m/s)  from  those  with  very  slow  (<2.5  m/s)  motion. 
An  attempt  along  this  line  was  made  and  alluded  to  in  Chapter  3. 

However,  the  method  of  separating  these  two  classes  of  cyclones  was  poor 
and  the  results  were  discouraging.  Better  results  could  probably  be 
achieved  if.  as  mentioned  in  Chapter  5,  a  larger  threshhold  value  of 


positive  or  negative  S00  mb  height  pattern  Indicator  was  used  and  the 
resulting  fast  and  slow  equations  were  used  on  only  those  cyclones  that 
exhibited  fast  or  slow  motion  based  on,  for  example,  their  past  12-hour 
motion.  Instead,  in  the  aborted  method  attempted  here,  all  cyclones 
used  either  the  fast  or  slow  equations  and  this  is  probably  where  the 
scheme  failed.  The  only  problem  with  the  suggested  substratification  is 
the  small  sample  size  that  might  result  for  the  very  fast  and  very  slow 
subsets. 

Other  improvements  may  be  possible  by  changing  the  grid  for 
computing  the  geostrophic  winds  from  the  geopotential  height  gradients. 
This  model  used  a  diamond  grid  with  grid  points  north,  south,  east,  and 
west  5°  away  from  the  point  where  the  steering  wind  was  being  computed. 
Patterns  other  than  a  diamond  or  distances  greater  than  3°  could  be 
tested.  The  use  of  actual  observed  winds  rather  than  height  gradients 
should  be  tested  although  the  use  of  wind  field  data  did  not  seem  to 
produce  any  slgnf leant  impact  on  forecast  errors  when  oompared  to  the 
use  of  height  field  data  for  Southern  Hemisphere  statistical-synoptic 
schemes  (Keenan,  personal  communication). 

In  any  event,  CSU84  -  as  developed  and  detailed  in  this  paper  - 
appears  to  offer  very  encouraging  verification  results  and  appears  to 
produce  significant  improvement  in  forecast  skill.  The  amount  of 
reduction  in  forecast  position  errors  suggested  by  these  results  should, 
hopefully,  hasten  the  testing  of  the  model  in  an  operational 


environment. 
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APPENDIX  A 


List  Of  All  Objective  Techniques  Currently  In  Operational  Use 

itic 

HURRAN  -  HURRicane  AN aloe  system  uses  analog  cyclones  in  a  time 

and  space  window  similar  to  the  current  cyclone  to  generate 
forecast  track.  (Hope  and  Neumann.  1970) 

CLIPER  -  CLImatoloev  and  insistence  system  is  a  purely 

statistical  technique  that  uses  climatology  and  persistence 
information  as  predictors  to  generate  regression  equations 
using  essentially  the  same  predictors  used  in  the  analog 
sense  by  HURRAN.  (Neumann.  1972) 

NHC67  -  statistical-synoptic  system  that  uses  linear  combinations 
of  observed  height  and  height  change  fields  at  the  1000. 
700.  and  500  mb  surfaces  as  predictors  for  zonal  and 
meridional  components  of  cyclone  motion.  Persistence  is 
also  used  for  the  early  forecast  periods.  (Miller,  e£  al. , 
1968) 

NHC72  -  combines  output  from  HURRAN.  CLIPER  and  NHC67  into  one 
forecast  using  a  modified  stepwise  screening  technique. 
(Neumann.  .§£  .al. ,  1972) 

NHC73  -  combines  output  from  CLIPER,  a  steering  subset  using 

height  gradients  and  thicknesses  and  a  synoptic  subset  that 
introduces  24-,  36-,  and  48-hour  500  mb  prognostic  height 
data  into  one  forecast  using  an  optimized  geographical 
stratification  system.  (Neumann  and  Lawrence,  1973) 

SANBAR  -  SANders  BARotroplc  model  is  a  filtered  barotropic  model 

using  observed  1000  to  100  mb  pressure  weighted  winds.  The 
wind  field  near  the  cyclone  is  modified  to  conform  to  the 
observed  cyclone  motion.  (Sanders  and  Burpee,  1968) 

MFM  -  Moveable  £ine  Mesh  barocllnic  model  has  10  vertical  levels 
and  60  km  horizontal  grid  spacing.  (Hovermale  and  Livezey, 
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TYphoon  ANalog  system  uses  analog  cyclones  with  similar 
time  and  space  characteristics  to  generate  24,  48,  and  72 
hour  forecast  positions.  (Jarrell  and  Somerville,  1970) 

CYcIods  500  is  an  updated  HATTRACK/MOHATT  steering  program 
that  provides  geostrophic  steering  forecasts  at  the  500  mb 
level  using  prognostic  fields  and  a  12-hour  persistence 
bias.  (Renard,  1968;  Renard,  et  al . .  1973) 

CLIMatology  aid  that  provides  forecast  positions  based  upon 
the  month  and  position  of  the  current  cyclone. 

Half  persistence  And  Climatology  provides  24-  and  48-hour 
forecast  positions  that  are  the  mid-points  of  straight  lines 
connecting  the  24-  and  48-hour  forecast  points  on  the  XTRP 
and  CLIM  tracks. 


EXTRaPolation  forecast  that  extropolates  a  straight  line 
from  the  12-hour  old  position  to  the  current  position. 

filended  persistence  &nd  Climatology  provides  24-,  48-,  and 
72-hour  forecast  positions  with  less  weight  being  given  to 
persistence  at  each  succeeding  forecast  interval. 

Hested  Tropical  Cyclone  Hodel  is  a  primitive  equation  model 
with  3  vertical  levels  and  a  41  km  grid  centered  on  the 
cyclone  that  moves  within  a  coarser  grid.  (Harrison,  1981) 

Cne-way  tropical  Cyclone  Model  is  a  three-layer  primitive 
equation  model  with  a  205  km  horizontal  grid.  Navy 
hemispheric  prog  fields  are  used  at  12-hour  intervals  to 
update  the  model's  boundaries.  (Hodur  and  Burk,  1978) 
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APPENDIX  B 

Correlation  coefficients  fields  for  the  Atlantic  on-the-ridge  and 
south-of -ridge  categories. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  on-the-ridge  fast  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(y)  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  on-the-ridge  slow  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(p)  and  the  stars  Indicate  the  gridpoints  at  which  predictors  were  selected. 


Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  south-of -ridge  fast  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(§)  and  the  stars  Indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  south-of-ridge  slow  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  luotion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(4)  and  the  stars  Indicate  the  gridpoints  at  which  predictors  were  selected. 


Correlation  coefficient  fields  (x  .01)  for  the  Atlantic  south-of-ridge  slow  V  notion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(6)  and  the  stars  indicate  the  grldpoints  at  which  predictors  were  selected. 
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Correlation  coefficients  fields  for  the  NW  Pacific  on-the-ridge  and 
south -of-ridge  categories. 
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Fig.  C.l.  Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  on-the-ridge  fast  U  motion.  The 

four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(4)  and  the  stars  Indicate  the  gridpolnts  at  which  predictors  were  selected. 


Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  on-the-ridge  fast  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (e)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(f )  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 


Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  on-the-ridge  slow  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
($)  and  the  stars  Indicate  the  gridpolnts  at  which  predictors  were  selected. 


to 
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change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure. 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
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Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  south-of-ridge  fast  U  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(6)  and  the  stars  indicate  the  grldpoints  at  which  predictors  were  selected. 
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Fig.  C.6.  Correlation  coefficient  fields  (x  .01)  for  the  NW  Pacific  south-of-ridge  fast  V  motion.  The 
four  panels  present  the  correlation  between  cyclone  motion  and  (a)  500  mb  height;  (b)  24-hour 
change  of  the  500  mb  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure. 
Northward  and  eastward  cyclone  motions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(£)  and  the  stars  indicate  the  gridpoints  at  which  predictors  were  selected. 
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Correlation  coefficient  fields  <x  .01)  for  the  NV  Pacific  soutb-of-ridge  slow  V  eotlon.  The 
four  panels  present  the  correlation  between  cyclone  notion  and  (a)  500  sb  height;  (b)  24-hour 
change  of  the  500  ab  height;  (c)  sea  level  pressure;  and  (d)  24-hour  change  of  sea  level  pressure, 
Northward  and  eastward  cyclone  notions  are  positive.  The  position  of  the  cyclone  is  denoted  by 
(4)  and  the  stars  indicate  the  gridpolnts  at  which  predictors  were  selected. 


Fig.  o.|.  Correlation  ooeffleleot  fields  Is  .01)  for  tho  forth  Indian  Oonb  on-the-rldge  0  notion.  The 
fwr  panel  a  praeent  the  oorr  elation  hate  a  on  cyclone  notion  and  (a)  900  nb  height;  (b)  24-hour 
change  of  the  300  nb  height;  {el  see  level  preeoure;  end  (d)  24-hour  change  of  oee  level  preesuo 
northward  and  eeatward  eyclooe  notlooa  are  poaltlve.  The  position  of  the  cyclone  la  denoted  by 
{£)  and  the  otara  indicate  the  grldpointa  at  which  predictors  were  selected. 
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iO  ab  streamline  analysis  for  1200  GMT  on  October  11.  Typhoon  Wynne  is  located  at 
approximately  300  n  mi.  east  of  Taiwan  at  this  time. 


(4)  24-hour  pressure  change  20°  south  20°  west  of  cyclone  ■ 
(2009-2020)  X  -8.2569  -  8.2S69 

(5)  pressure  20°  south  of  oyolone  •  2006  x  3.4909  -  3322.8434 
(4)  24-hour  pressure  change  20°  north  3°  west  of  oyolone  ■ 

(1027-1024)  x  -1.4298  -  -4.2894 

(7)  24-hour  height  ohsnge  13°  north  23°  west  of  oyolone  ■ 
(3742-3734)  x  -.2148  -  2.3776 

(8)  24-hour  pressure  ohsnge  3°  north  13°  east  of  oyolone  ■ 
(1016-1017)  x  4.6907  -  -4.6907 

(9)  pressure  3°  north  3°  west  of  oyolone  -  1013  x  -3.9711  - 
-4030.6663 

(10)  height  10°  south  40°  west  of  oyolone  -  3844  x  .2906  ■ 
1698.2664 

(11)  past  24-hour  east-west  notion  “  1111.949  x  oos  (21.3  +23.8)/2) 
(127.2-130. 8) /24  -  -133.9284  x  -.1238  -  19.0363 

(12)  24-hour  height  ohsnge  30°  north  3°  east  of  oyolone  - 
(3363-3329)  x  -.1294  -  -4.3996 

(13)  pressure  20°  east  of  oyolone  »  1011  x  4.7028  -  4634.3308 

(14)  24-hour  pressure  change  13°  north  40°  west  of  cyclone  - 
(1018-1018)  x  1.6731  -  0 

(13)  24-hour  height  change  13°  north  20°  east  of  oyolone  - 
(3763-3777)  x  -.1196  -  1.6744 

(16)  24-hour  height  ohsnge  33°  north  30°  east  of  oyolone  - 
(3186-3161)  x  -.0761  -  -1.9023 


(17)  24-hour  change  of  tho  north-south  ooaponont  of  tho  200  ah  steering 
over  ths  oyolons  ■ 

([ (12407+12489+12499*12489*12479)/ 3  -  ( 12372  *12454*124<8*12464*12460)/ 
5)  /  (sin (23.9)  z  oos(23.8>)  -  I(12423*12487*1230S*12491*12480)/S  - 

( 12376+12462+1247 6+12470+12469)/ 3]  /  (sln(23.t>  Z  OOS(23.t)))  z  .0499 
-  (71.9423  -  73.1236)  z  .0491  -  .2494 

(18)  24-hour  bslfht  change  33°  north  30°  west  of  oyolons  • 

(3471-3488)  z  -.0844  •  -1.4348 

(19)  24-hour  Chang*  of  ths  sast-vsst  ooaponont  of  ths  200  ah  stssrlng 
10°  north  of  cyclone  - 

[ (12483-12208) /aln( 33. 8)]  -  [ (12479-12187 )/sin(33. 8)]  z  -.0333  - 
1.0787 

(20)  24-hour  prssaurs  ehangs  10°  south  10°  sast  of  cyclons  - 
(1011-1011)  z  -3.98  «  0 

(21)  Intsrespt  -  -4704.3091 

Sunning  all  of  ths  contributions  (103.3249)  and  adding  it  to  ths 
initial  latltuds  results  in  a  24-hour  forsoast  latituds  of 

((103.3249  z  24)  /  1111.949]  *  23.8  -  24.0 

Applying  ths  rogrssslon  equation  for  ths  sast-wsst  dlsplaeensnt  for 
ths  on-ths-rldgs  oatsgory.  ths  first  prsdlotor  is  tbs  300  ab  height  10° 
north  and  3°  vast  of  tbs  oyolons  center.  Proa  Table  £.1.  this  height  is 
3842  a.  Multiplying  this  height  by  its  oo efficient  -.464  results  in  s 
contribution  of  -2710.688.  The  ssoond  prsdlotor  is  ths  past  12- hour 
east-west  notion.  Ths  initial  oyolons  position  is  23.8°N  127. 2°E.  Ths 
past  12-hour  position  is  22.6°N  128. 8°E.  Ths  past  12-hour  east-west 
dlsplaosasnt  is  oonputed  as: 


I  -  1111. Mi  x  oosldl. <  ♦  2S.t)/2]  x  (1272-120. 0)/12 
-  -124.2 

Multiplying  fcblx  dlaplaoeaeat  by  lta  ooof flcloot  .202  rooultx  lo  o 
contribution  of  *41.142.  Slollxrly.  prodlotorx  (2>  through  (20)  rooultx 
lo: 

(2)  prooauro  10°  xouth  of  oyoloao  -  1004  x  4.497  ■  4222.902 

(4)  height  10°  north  12°  oxat  of  oyelooo  ■  2071  x  -.241  ■ 

-1222.221 

(2)  24-bour  height  change  2°  north  l0°  weot  of  ay  clone  - 
(2041-2074)  x  -.420  -  0.144 

(4)  24-bour  proa  aura  ah cage  22°  north  22°  Meat  of  eyelone  • 

(1029-1022)  x  1.22  -  -4.29 

(7)  prooauro  2°  aouth  22°  oaat  of  oyoloao  ■  1011  x  2.021  •  2044.241 

(0)  24-hour  prooauro  change  10°  aouth  20°  Moot  of  oyoloao  - 
(1007-1000  x  0.772  -  -0.772 

(9)  height  40°  oaat  of  oyoloao  •  2002  x  .249  •  2022.010 

(10)  24-hour  height  ohango  10°  aouth  10°  Moat  of  oyoloao  • 

(2029-2070)  x  -.092  -  9.012 

(11)  24-hour  height  ohango  20°  north  12°  Moat  of  oyoloao  • 

(2474-2211)  X  .114  -  -4.210 

(12)  24-hour  prooauro  ohango  10°  north  20°  oaat  of  ay clone  • 

(1017-1022)  x  -2.920  -  14.44 

(12)  prooauro  20°  north  22°  oaat  of  oyoloao  ■  1000  x  1.92  •  1920 

(14)  paat  24-hour  north-aoutb  notion  -  (1111.949  x  (22.0-21.2)/24)  x  .102 
-  19.2 

(12)  oaat-Moat  oooponent  of  the  200  ab  at  earing  over  the  oy alone  ■ 

[(12 440+12 442 +12 471 +12 47 9)/ 2  -  (12272+12292+l2411+12412+12407)/2)  / 
aln(22.9)  x  .042  -  10.02 

(14)  24-hour  prooauro  ohango  22°  north  22°  oaat  of  oy alone  ■ 

(997-904)  x  -1.222  -  -14.422 

(17)  height  22°  Moot  of  oyoloao  -  2049  x  -.219  -  -1042.021 

(10)  height  2°  aouth  10°  Meat  of  eyelone  •  2044  x  -.249  -  2144.224 


(It)  u  I wr  boi^t  «fe4B|«  IS0  aortt  IS°  tut  of  «relOM  • 

(2402-2202)  s  -.022  •  -.717 

(20)  24-boor  b«l«bt  IS4  oortb  21°  v«t  of  cpeloao  • 

( 5444-54* t)  >  -.221  •  .204 

(21)  Xotoroopt  •  -7002.245 

I— if  oil  tfei  ooBtrlbotlooo  (-20.222)  ood  oddlog  It  to  tbo 
initial  lnogifdt  rooulta  lo  •  24-boor  foroooot  loogltodo  of 

(-20.222  M  24)  /  1111. Mt  l  ooo((22.f  ♦  20.0)/2)  ♦  127.2  •  122.0 

la  o  util  oodo.  grltfpolat  dot*  Mould  oot  bo  oroiloblo  to  tbo 
feroeootor.  Xootood.  bo  soot  ooo  tbo  atoUoblo  oboonrotlooo  ood 
ootlaoto  bolgbta  ood  oroooaroo  ot  tbo  roqolrad  grldpo&oto.  To  toot  tbo 
opnojtlTlty  of  tbo  roproooloa  oqoatlooo.  oil  of  tbo  roqoirod  boltfito 
moto  oltbor  loorooood  or  doorooood  by  20  o  froo  tbooo  tbot  appear  lo 
Toblo  1.1.  Similarly.  oil  oT  tbo  prooouroo  ooro  oltbor  loorooood  or 
doorooood  by  2  ob  froo  tbooo  tbot  oppoor  la  Toblo  K.2.  ■  oltbor  tbo 
toodooeloo  oor  tbo  otoorloj  prodlotoro  ooro  oltorod  olaoo  It  ooo  folt 
tbot  24-boor  obootoo  ood  baMbt  grodlooto  oro  ooolor  to  ootlaoto  tboo 
tbo  oboolato  trolooo. 

for  tbo  oooo  of  loorooolog  oil  or  tbo  botpbto  ood  proooaroo.  tbo 
foroooot  lotltodo  woo  bolf  o  do po  fartbor  oor  to  ood  tbo  foroooot 
loatltodo  ooorly  o  dotroo  fortbor  ooot.  l.o. .  2(.2°b  ood  122. 0°t.  for 
tbo  oooo  of  door ooo log  oil  of  tbo  boifbto  ood  proooaroo.  tbo  foroooot 
lotltodo  woo  bolf  o  dtp>M  fortbor  ooatb  ood  tbo  foroooot  loagltado  woo 
ooo  dagroo  fartbor  woot,  l.o..  2J.5°i  tad  122. 0°K.  Tboo.  tbo  worot  oooo 
oltootloo  of  all  tbo  bUptto  ood  proooaroo  balog  too  blpk  or  too  low 
weald  ’odd'  opproxlootaly  20  o  al  to  tbo  foroooot  poo It loo.  lotvalot 
tbot  20  a  for  tbo  200  ab  bodfbt  ood  2  ob  for  tbo  oorfooo  prooooro  oro  o 


r«UiU«  r«i>  af  valNi  UMt  a  aaalyaia  mil  vary  ma  tuv 

crlApalat  iMliala.  •  M  »al  Hffar—  U  tka  far— aat  pmUIn  m«U 
raadt. 

Two  palata  akaal*  ka  aapkaataad.  Tka  flrat  to  tkat  ralatlva  to  tka 


alayia  aotfal  —an  aa  tkla  la  tka  oyyartaalty  far  tka  faraaaatar  to  layot 
tka  data  vita  aakjaatlva  taalltf  awtrai.  If  tka  apaaytla  valaw  at 
yridpolata  raqalrad  ay  tka  aodal  ayyaar  faaatlaaakla.  tka  faraaaatar 
aaald  aaally  aailfy  tkaa  ta  sera  raaUatla  valaw.  ttltk  tka  aid  af  a 
daak  toy  aaloalatar.  ka/aka  aoall  raroa  tka  wdal  aovoral  tl—  v&tk 
dlffaraat  pradlator  valaw  ta  aw  tka  affaot  of  tka  *wtw.  Paul  tka 
aatpat  fraa  tka  ylokal  Pt  aadol  lapravw.  aakjaatlva  aodlfloatlw  of  tka 
feraawt  fialda  aaak  w  oat  1  load  kora  will  prakakly  ka  aaaaaaary. 

A  atatlatlaal  Ayawlaal  traok  pradlatlw  aodal  la  oat  aarraatly 


avallakla  for 


la  tka  Mf  Paalfla/Swtk  Gkiw  Aw 


tM 

i.  la  mrnm  fev  tta  tarmdaUw  af  tM  Mtel  la  Clatter  I.  tM 
Ual  far  Hpl/lMtt  lapnMMt  art  at  a  with  tM  im  af  ttla  aaial. 
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